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ABSTRACT 

The EFAR project is designed to measure the properties and peculiar motions of 
iy) | early-type galaxies in two distant regions. Here we describe the maximum likelihood 

algorithm we developed to investigate the correlations between the parameters of the 
| EFAR database. One-, two-, and three-dimensional gaussian models are constructed 

to determine the mean value and intrinsic spread of the parameters, and the slopes and 
intrinsic parallel and orthogonal spread of the Mg 2 -Mg6', Mg 2 -cr, Mg&'-cx relations, 
and the Fundamental Plane. In the latter case, the cluster peculiar velocities are 
also determined. We show that this method is superior to "canonical" approaches of 
least-squares type, which give biased slopes and biased peculiar velocities. We test 



the algorithm with Monte Carlo simulations of mock EFAR catalogues and derive 
•£3 ■ the systematic and random errors on the estimated parameters. We find that random 

^ I errors are always dominant. We estimate the influence of systematic errors due to 

the way clusters were selected and the hard limits and uncertainties in the selection 
function parameters for the galaxies. We explore the influence of uniform distributions 
1 in the Fundamental Plane parameters and the errors. We conclude that the mean 

' peculiar motions of the EFAR clusters can be determined reliably In particular, the 

placement of the two EFAR sample regions relative to the Lauer & Postman dipole 
allows us to strongly constrain the amplitude of the bulk motion in this direction. We 
justify a posteriori the use of a gaussian modeling for the galaxy distribution in the 
Fundamental Plane space, by showing that the mean likelihood of the EFAR sample 
is obtained in 10 to 30% of our simulations. We derive the analytical solution for the 
maximum likelihood gaussian problem in N dimensions in the presence of small errors. 

Key words: galaxies: clusters - galaxies: elliptical and lenticular, cD - galaxies: 
fundamental parameters - - galaxies: distances and redshift - cosmology: large scale 
structure of Universe 



1 INTRODUCTION 

The EFAR collaboration (Wegner et al. 1996; Paper I) has collected photometric (Saglia et al. 1997a, 1997b; Paper III 
and IV) and spectroscopic (Wegner et al. 1999; Paper II) data for galaxies in clusters with the primary goal of using the 
tight correlations between the global properties of early-types (the Fundamental Plane, hereafter FP, and the D„-a relation, 
Djorgovski & Davis 1987, Dressier et al. 1987) to measure the peculiar motions and the mass distribution on large scales. 
Colless et al. (2000, Paper VII) discuss the results of these investigations and the interpretation of the peculiar motions in the 
context of cosmological models. This paper presents the methodology adopted to derive these results and the Monte Carlo 
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simulations performed to test them and estimate the uncertainties and residual systematic biases. Colless et al. (1999, Paper 
V) discuss the implications of the EFAR Mg-<r relation. 

An abundant literature is dedicated to the problem of how to best determine regression lines from incomplete datasets 
subject to errors and explicit selection criteria. Regression algorithms involving two variables are discussed by Isobe et al. 
(1990) and Feigelson & Babu (1992), who provide least-squares fits and their uncertainties. Akritas & Bershady (1996) address 
the problem of linear regressions with errors and intrinsic scatter, by considering the variance of the variables involved. La 
Barbera, Busarello and Capaccioli (2000) extend this work to the analysis of the FP equation. The canonical methods used 
to find the coefficients of the FP (or of correlations between two variables) are of the least-squares type: one minimizes the 
orthogonal (absolute or squared) residuals from the plane (j0rgensen et al. 1996, J0rgensen 1997, Pahre et al. 1997, Scodeggio 
et al. 1997) or the (squared) residuals from one of the variables (Hudson et al. 1997). Principal component analysis has also 
been used (Carvalho & Djorgovski 1992). Ample discussion of problems related to selection is given by Willick (1994) and 
Teerikorpi (1997), in the context of the determination of peculiar velocity fields. However, none of these approaches is able to 
deal effectively with the multiple problems that the EFAR and similar datasets pose: (i) the factor two in redshift spanned, 
(ii) the strong selection effects, (iii) the non-negligible and widely varying measurement errors, (iv) the presence of defined 
selection criteria, (v) the intrinsic scatter of the relation. In this paper we demonstrate that "classical" methods fail when 
confronted with the EFAR dataset, and we quantify the biases they produce. We construct an algorithm based on Maximum 
Likelihood (ML) gaussian modelling and test it with Monte Carlo simulations. We demonstrate that it solves the problems 
posed by the dataset, producing nearly unbiased estimates of the parameters involved. 

The paper is organized as follows. ^ discusses critically the results of linear regression analysis and presents the ML 
algorithm for the one-, two- and three-dimensional cases. ^ describes the method to generate mock catalogues of the EFAR 
database, the tests of the "canonical" least-square methods and of the ML gaussian algorithm, and the estimates of the 
errors on the derived parameters. Simulations backing the results presented in Papers II, V and VII are discussed for the 
Mg2-Mg&', Mgfe'-cr, Mg2-cr relations, the Fundamental Plane and the cluster peculiar velocities. Conclusions are drawn in 
The Appendices give the analytical solution of the N-dimensional gaussian maximum likelihood problem in the presence of 
small errors and some results used in the paper. 



2 THE MAXIMUM LIKELIHOOD GAUSSIAN ALGORITHM 

After reviewing the results of the standard least-squares analysis in 



2.2 



2.3 



2.1 



we describe the properties of the EFAR database in 
to argue that a ML approach is needed to st udy the sa mple . We set th e eq uations of the problem for the general case in 
and address the simplified cases when one (§2.4), two (§2.5) or three (§|2.6|) variables are considered. 



2.1 Linear models 

Before discussing the specific problems related to the EFAR database and the solutions we have developed to address them, 
it is worth considering the more general question of how to model the correlations existing between the structural parameters 
of galaxies. As usual in astrophysics, logarithmic quantities are considered in order to avoid scaling problems. In the simplest 
approach, one focuses on pairs of datapoints {xi, yi}, looking for a linear relation y = ax + b between them. Table 1 of Isobe et 
al. (1990) (see also Eqs. below) summarizes the estimates of the slopes a and zeropoints 6, and their variances, derived using 
five methods (Y-X and X-Y regressions, bisector, reduced major axis and orthogonal regression), when measurement errors 
are negligible. Briefly, we recall that the Y-X regression derives the slope ayx and zero point byx of the line y = ayxx + byx 
by minimizing the quantity: 

XYX = — dYXXi — byx) 2 , (1) 

and the X-Y regression by minimizing the analogous Xxy- The bisector solution gives the slope and zeropoint of the line 
passing in between the Y-X and X-Y lines. The reduced major axis slope is the geometric mean of the Y-X and X-Y slopes. 
The orthogonal regression derives ao and bo by minimizing the orthogonal residuals from the line: 

2 _ v (Vi -aox l -bp) 2 
1 + "o 

If the differences in slopes and zeropoints between the methods are smaller than the expected variance, a well-defined answer 
is found. Unfortunately, this is rarely the case in astronomy. Therefore, Isobe et al. (1990) recommend astronomers use the 
Y-X regression when it is clear that Y is the variable to be "predicted" , and the bisector method when the functional relation 
between the variables has to be investigated. They discourage the use of the orthogonal regression, since its slope has greater 
dispersion. The problem is complicated when errors (on one or both variables) are present. In particular, Feigelson and Babu 
(1992) conclude that "models incorporating both measurement errors and intrinsic scatter are complex and not yet fully 
developed". Akritas and Bershady (1996) improve on this aspect, presenting an updated version of Table 1 of Isobe et al. 
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(1990), where corrections for measurement errors are taken into account. Isobe et al. (1990), Feigelson and Babu (1992) and 
Akritas and Bershady (1996) stress that "there is no such thing as true slope", and that the regression method should be 
chosen according to the problem to be solved. 

These considerations become more transparent when the underlying bivariate probability distributions P(x, y) are exam- 
ined. Assuming normality, P(x,y) can always be written as: 



P(x,y) 



■. exp 



2(1 -P 2 ) 



(x-x) 2 2p(x-x)(y-y) (y - y) 2 



(3) 



where x and y are the mean values of x and y, a x and a y the standard deviation (rms scatter), and p the correlation coefficient 
defined as o xy = pa x a y . Calling x — x — x and y = y — y , Eq. ^| can be cast as: 

\A\^ 2 



P(x,y) 



2tt 



■ exp j-- [A xx x 2 + 2A xy xy + A yy y 2 ] j 



where A = V \ |A| = A xx A y 



A z = 

ll - Xy 



(1 — p 2 )a 2 a 2 ] 1 and V is the covariance matrix: 



V = 



P0 x 0y 



pO X Oy 
2 



A = 



Vol 



1 - P 2 V -p/( a * a y) 



-p/(a x a y ) 
l/al 



Eqs. ^ and ^ are equivalent to the three following formulae: 



P(x,y) 

P(x,y) 

P(x,y) 
where: 



2lXO X 



27TCT, 



■ exp 



■ exp 



x 

2aJ 



2o 2 



2%OyyfT 



: exp 



i [y 



^k<Tx yl - p : 



: exp 



27TCT1 



■ exp 



(y - a x) 2 
~2{l + al)al 



2-KU2 



■ exp 



2 (l-p 2 )a 2 

2 (l~p 2 )a x 
(a y + x) 2 



2(l + a2)»j2 



ol + 



2 
""I 



2pcr a; crj / 
i(aj + a 2 y + ^JJpl 



\J (pi - a 2 ) 2 + 4:p 2 a 2 a 2 



o 2 Y + 4f?oio$) 



(4) 
(5) 

(6) 

(7) 

(8) 

(9) 
(10) 
(11) 



°l = \(ol + al - ^(a 2 -a 2 Y + ^ala 2 ). 

Eq. ^ shows that a gaussian bivariate probability distribution can be generated extracting first the x variable around its 
mean with rms a x , then the residual £ = y — ^-x around the Y — X line with slope equal to the regression of Y on X and 
rms yjl — p 2 a y . Similarly, Eq. fi] shows that the same distribution can be generated first extracting the y variable and then 
the residual f = x — ^-y around the Y — X line with slope equal to the regression of X on Y. Finally, Eq. ^ generates the 
distribution extracting the residuals around the eigenvectors of the covariance matrix. The direction of the eigenvector with 
the largest eigenvalue is the slope of the orthogonal regression. The zero points b of the relevant linear relations of slope a are 
obtained from the mean values as: 

b = y~ax. (12) 

Do we therefore conclude, as Akritas and Bershady (1996), that "there is no such thing as true slope"? In fact, the key 
of the problem is the covariance matrix, and therefore the natural slope is the direction of the principal axes (given by Eq. 
^|) with the largest eigenvalue (given by Eq. |lo[ Eq. [n] is the second eigenvalue of the covariance matrix) . As a consequence, 
the orthogonal regression determines the "true slope", not the bisector, which Isobe et al. (1990) prefer for its lower variance. 
However, what about distance determination problems, where Eq. |^ should be preferred according to Akritas and Bershady, 
if y is the distance-dependent quantity? Indeed, if one seeks for the distance shift 8 of the datapoint (j/; + 5, Xi), the most 
probable value when Eq. ^ is given is: 



•(Vi - —Xi), 
a x 



(13) 



which involves the Y — X regression coefficient V xy /V xx . 

Further complications arise when the datapoints {xi,yi} are affected by errors comparable to the rms spreads a x and 
a y , or when only a subset of the possible data volume is available due to incompleteness. These affect the estimation of the 
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covariance matrix V and the mean values of the variables. For example, we estimate in Appendix ^ that the orthogonal 
regression solution ignoring the errors underestimates the true slope if the error in the X direction is larger than the one in 
Y. These effects might be included in the regression approach by careful treatement (see references given in Feigelson and 
Babu 1992), but it is clear from the discussion above that a maximum likelihood approach is superior, allowing the natural 
and simultaneous solution of all these difficulties. 



2.2 The EFAR database 

The EFAR database comprises photometric and spectroscopic information for a set of cluster galaxies spanning a factor of 
two in redshift, with non-negligible measurement errors and selection effects. 

The size of the errors can be quantified comparing the rms spread of the variables with the mean of the errors. In this 
sense photometric errors (« 0.02 mag, see Paper III and IV) are small in the EFAR database. Unfortunately, this is not the 
case for the spectroscopic data (see Paper II). For the case of the velocity dispersions a, we find that the mean value (25 
km/s) of the measurement errors on a is 36% of the rms spread of the values (71 km/s). For the Mgb index, the mean value 
(0.39 A) of the measurement errors is 48% of the rms spread of the values (0.83 A). In addition, we cannot measure central 
velocity dispersions smaller than 100 km/s for lack of spectral resolution. 

Selection effects are also severe. The galaxy sample is more than 50% complete only for galaxies with Dw > 20 kpc (see 
Paper I), with 46% of the galaxies having selection probabilities larger than 50%. The EFAR galaxy selection S is described 
in Paper I and is a function of the Di : w diameter measured in arcsec for galaxy i: 

5<j(logA,w) = 0.5 1 1+ erf 

where j is the cluster index, log Dyyj is the midpoint and 5w,j the width of the cutoff in the selection function. The diameter 
Di t w correlates with the Di <n diameter of the sample as log Di, n = 0.8 * log Di : w + 0.26 with 0.09 dex scatter in Di <n (see 
Paper III). Therefore we compute the selection probabilities Si,j from the Di t wn diameters derived from the very accurate 
Di, n instead of using the values of Di,w given in Paper I. In the following we shall indicate with Si the selection probability 
for the galaxy i, dropping the cluster index, and with Wi = 1/Si the selection weight. 

Due to the large spread in redshift, the EFAR clusters have been sampled down to different limiting D^(kpc). However, 
95% of the early-types of the sample have Dwn > 12.6 kpc, so that nearly unbiased estimates of the cluster peculiar velocities 
can be obtained just by dropping from the sample objects with \ogDwn < d cu t = logDwcut- 

To our knowledge, no algorithm based on linear regressions has been developed to model all these features at the same 
time. In the following we describe a parametric algorithm based on gaussian modelling able with to deal with (i) a spread of 
selection weights, (ii) sizeable measurement errors with a large spread, (iii) possible explicit selection limits. 



log D itW - log D w j 



OW,j 



(14) 



2.3 The general case 

In the general case, we have N (logarithmic) data xl = (xi,i, ■ ■ ■ , JEj.jv) for each galaxy i. We indicate with Ei the related 
N x N error matrix. If the errors are uncorrelated, Ei is a diagonal matrix with Eij,k = 5j,fc<xf where 8j,k = 1 if j = k and 
otherwise. Each galaxy has been selected according to the value of a selection diameter, which possibly could be expressed as 
a linear combination of a subsample of Xi t N, giving a selection weight w% = 1/ Si. One variable Xi^\ may be distance dependent. 
We want to estimate the vector of the mean values of the data x — (xi, ■ ■ ■ ,~xn), the covariance matrix V and, possibly, the 
vector S — (Si, ■ ■ ■ , 8m) of the M cluster peculiar velocity shifts. In this case we need an additional constraint to fix xi, such 
as ^2 ■ 8j — 0. Other options are discussed in Paper VII. The direction of the eigenvector with the smallest eigenvalue of the 
matrix V defines the minimum variance hyperplane describing the data distribution, or the slope of the linear correlation in 
the two dimensional case. We solve the problem assuming that both the distributions of the galaxy parameters and of their 
errors are normal, with covariance matrix V and the error matrices Ei respectively. The first assumption can be verified a 
posteriori, estimating the likelihood of the best-fitting gaussian model (see Figs. ^, [l3| and §3.5.3). The second is based on the 



error analysis performed in Papers II and III for the EFAR database. Once the error convolution is performed, the probability 
density of the vector Xi is: 

= ,„ w/21,/ 1 , tFTUWI ex P f-o^V + Ei)- 1 ^] §{Aii - x cut ) (15) 



(2tt) n / 2 \V + E t \y^fi 



2 



where Xi = (xi,\ —xi + 5j, Xi,2 — X2, ■ ■ ■ , x%,n — xn)- The function 9(y) — Y[ k S(yk), where 9(y) — 1 if y > and otherwise, 
takes into account that parts of the parameter space might not be accessible because of selection effects or explicit cuts. For 
simplicity we assume that these cuts are applied to linear combinations of the variables, described by the appropriate matrix 
A. The normalization factor /; is such that j Pd N x — 1. Following Eadie et al. (1971), we write the likelihood of the observed 
sample as: 
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c^llp&r, (i6) 

i 

counting Wi = 1/Si times galaxies with selection probability Si We determine x, V and, possibly, S by minimizing — ln£. In 
the following sections we describe how this general scheme is implemented in the different cases of interest. In all cases the 
minimization is performed numerically, using the simplex algorithm (Press et al. 1986). Appendix ^ discusses the analytical 
solution of Eq. hq for small errors, no explicit cuts and no peculiar velocities. 



2.4 The 1-dimensional Case 

The first simple application of the scheme described above is the determination of the mean value p and the (intrinsic) rms a 
of a set of data a, with errors <jj, selection weights w, and subject to the censoring x > x cut . Mean values and intrinsic spread 
of the EFAR photometric (logRe, (SB e ), \ogD n , (B — R)) or spectroscopic (logo - , Mg6, Mg2) parameters, or of the measured 
cluster peculiar velocities can be estimated in this way. 
The simplified version of Eq. |l(] reads now: 



C 



n 



fuy/2ir(cT? + a 2 ) 



exp 



M) 2 



6{xi 



Xcut ) 



(17) 



Errors are taken into account by convolving the intrinsic distribution with the appropriate gaussian error distribution. Selection 
effects are taken into account in Eq. [j^ by counting 1/ Si times galaxies with selection probability Si . The presence of explicit 
cuts Xcut is taken into account (for mild cases; see below) by the normalization correction: 



fu = 



0.5 | 1 - erf 



C, exp (- 



tlx 



(18) 



We expect our approach to work in the presence of mild cuts with (x cu t — p)/a << 0, so that the cut is well away from the 
mean value of the distribution (see §3.2). Appendix [X] gives the analytical solution of the problem for the limiting case of 
small errors. 



2.5 The 2-dimensional Case 

We turn now to the two-dimensional case, which allows us to investigate correlations between pairs of variables, such as 
Mg2-Mg6' examined in Paper II and the Mg&'-o and Mg2-o discussed in Paper V. Other relations that we shall not explicitly 
consider here are the FP - logo (where FP = logR e ~ b(SB e ) with b « 0.3), the (B - R) - logo-, (B - R) - Mg& and 
(B — R) — Mg2 relations, the "Kormendy" relation log R e — {SB e }, or the log D n -log a relation. Note that the EFAR database 
measures the photometric parameters in the R band. 
Eq. [l5| reads in this case: 

|AJ 1/2 1~t 

P(xi = (xtifXi*)) = tt - ? — exp(--f i (V + Ei) Xi)9(xiA - xicut)0(xi,2 - x 2c ut), (19) 

i7T/2i 2 

where the normalization factor f 2 i is 

fx = f^J^expi-^A^dxdy (2Q) 

= L(h,k,p). 

Here Aj = (V + Ei)^ 1 and L(h, k,p) is the bivariate probability integral (Abramovitz and Stegun 1971) with h = (xi cu t — 
xi)-y/Aj,22(l — p 2 ), k = (x2cut — X2)y/Ai,u(l — p 2 ) and p = — Ai,i2/^/Ai,nAi, 2 2. We set xi cu t = — oo, x 2c ut = — oo (and 
therefore f 2 i = 1) when considering the the Mg2-Mg6' relation, xi cu t = — oo, x 2c ut = logo cu t, when studying the Mg&'-o and 
Mg2-o relations, and x\cut = FP cut , x 2cu t = log a cu t, when examining the FP — a relation. In this case we derive FP cut using 
its relation with Dwn'- FP = 0.78 log Dwn — 6.14 (see Paper III). Finally, note that in every case we compute A, using the 
simplifying assumption that the error matrix is diagonal with diagonal terms given by the estimated total errors. For the 
spectroscopic data these involve two terms, the uncorrelated statistical errors and the correlated errors coming from run-to-run 
corrections (see Paper II). 

Appendix M gives the analytical solution for the limiting case of small errors, no peculiar velocities and no cuts. 
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2.6 The 3-dimensional Case 



The three-dimensional case allows us to study the EFAR Fundamental Plane. As is well known, early-type galaxies do not 
fill the three-dimensional space denned by the coordinates x = (x\ = log_R e ,a;2 = logcr,X3 = (SB e ) ), but rather occupy a 
narrow region around the Fundamental Plane defined by the equation: 



log R e — alog a — b(SB e ) = c. 



(21) 



Following §2.3, we determine the mean values log-Re, logo" and (SB e ), the covariance matrix V and the cluster peculiar 
velocity shifts Sj (subject to the constraint £\ 8j — 0, see §2.3) by maximizing the likelihood C: 



ln£ 



E 



T>CTcut,FP>FP C ut 



(log 14 



■ In P(xi) 



E 



a>lT c „t,FP>FP cl 



0.5^ (V + Ei)- 1 ^ +ln/ 3i + 1.51n(27r) +0.51n|F + E t 



where 



/3i 



PiOQog a - log a cut )6(FP - FP cut )d 3 x 



(22) 



is the fraction of galaxies with a > o cut and FP > FP cut . The integral (Eq. |22| ) is performed in Appendix The offset Sj 
between true mean galaxy size log_R e , and the mean galaxy size observed for cluster j, logR e — Sj, is related to the peculiar 
velocity of the cluster. In particular, the ratio of the true angular diameter distance of a cluster, Dj, to the angular diameter 
distance corresponding to its redshift, D(zj), is: 

Di 



dex(log R e ) 



(23) 



D(zj) dex(logi? e -5 

The relation between angular diameter distance and redshift (Weinberg f 972) is given by 

[Z> ~ H (l + z) 2 1 + qo z + VI + 2g 2 ' 

We assume _ffo=50kms _1 Mpc, go=0.5, and compute all redshifts and peculiar velocities in the CMB frame of reference. The 
peculiar velocity of the cluster, Vj, is then obtained as 

-zj - cz(Dj) 



(24) 



Vj ~ l + z(D 3 ) ' 

where z(Dj) is the redshift corresponding to the true distance Dj through the inverse of Eq. H 



(25) 



The error matrix Ei is given by: 



Ei = 



( 



\ 



8r\ 






2 





(l + a 2 )Sr 2 -SFP 2 





-l)SFP 2 + (l + a 2 )Sr 
a 2 (l + a 2 ) 



(26) 



+ 8ZP? 



where Sri is the error on log-R e ,i, SFPi is the error on the combined quantity FPi = logi? e ,i — a{SB e ,i), with a ~ 0.3, and 
SZPi is the photometric zero-point error (see Paper III and IV). Note that in principle b might be different from a. The 

quantity 8ui = \J ^° 1 ^ S ^Z(i+a* 2 ~) C " ^ s error on t ne effective surface brightness coming from the fitting. See Appendix 
^ for a derivation of these formulae. The error a s on log a combines the uncorrelated and correlated errors (see Paper II). 

It is worth noting that the eigenvectors v{, 1T2, V3 of the V matrix can be written to a close approximation as a function 
of the parameters a and b defining the Fundamental Plane of Eq. pit 



Vi = x± — aX2 — bx:j, 
V2 = Xi + xz/b, 



(27) 



V3 



= — Xi/b — (1 + b )x 2 /(ab) + x 3 , 



where vi measures the displacement from the FP, and V2 and V3 are two orthogonal vectors in the FP. The absence of a X2 
component in V2 is rema rkable . In Paper VII (Eq. 12) we fit for the coefficient of log a in the second principal component and 
find that it is small. In 



5.2 (Eq. BOh we show that it is compatible with zero within the statistical uncertainties. 



In the following we indicate the rms spread around the vectors vi , V2 and A3 with 01, cr% , cr 3 , respectively. 
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3 MONTE CARLO SIMULATIONS 



The algorithms described in the previous section have been extensively tested on mock catalogues of the EFAR database. 



After describing how these catalogues are generated (3.1), we proceed to test cases of increasing complexity. We start with the 
one-dimensional a distribution (3.2), as an example of the influence of hard cuts on the mean and rms of the sample. Then 
we consider the Mg2-Mg&' relation (|3.3|), where distance-independent quantities not subject to explicit cuts are involved. We 
turn to the Mg-cr relations (3.4), where no datapoints with central velocity dispersions smaller than the resolution limit can 
be present. Finally, we examine the three-dimensional case of the Fundamental Plane (3.5), where one distance-dependent 
quantity (log R e ) is correlated against two distance-independent quantities (log a and (SB e )), with distance-dependent selection 
limits and the presence of cluster peculiar velocities. 



3.1 Mock catalogues of the EFAR database 

We generate mock catalogues of the EFAR database using the gaussian distribution functions that are reconstructed from 
the data (see Paper II, Paper V, Paper VII and Case of Table 0). We specify the values of logi? e , log cr, (SB e ), Mg 2 , Mgb', 
the coefficients of the FP and the slopes of the Mg2-<r and Mgb'— a relations, and the dispersions of the distributions in the 
orthogonal directions. We generate one entry in the mock catalogue for each early-type galaxy present in the database. As a 
first step we extract the triplet (log-R e , logo - , (SB e )), using Eq. ^to specify the directions of the principal components of the 
distribution in the (log-R e , logo", (SB e )) space. We determine the corresponding value of the D n diameter using the relation 
(a second order approximation to the equation 4 of van Albada, Bertin and Stiavelli 1993): 

log D n = log 2R e - 0.289 A - 0.019 A 2 , (28) 

where A = (SBe)— 20 matches the surface brightness used to determine the EFAR D n diameters (see Paper III). Only galaxies 
with a dominant exponential component in the luminosity profile deviate strongly from Eq. (pj[). We simulate the effects of 
cluster peculiar velocities by adding cluster-dependent shifts to the extracted values of logJ? e and logD n . We convert \ogD n 
into log Dw using the relation log Dw = (logD n — 0.27)/0.8 (see Paper III) and add the 0.09 dex rms random scatter of 
the relation. We compute the corresponding selection probability, according to the selection parameters of the cluster being 
extracted and its (redshift) distance and test if the extracted datapoints should be in the sample. If this is the case, the values 
of Mg2 and Mgb' are extracted following the (two-dimensional) gaussian distribution computed at logo". Measurement errors 
are assigned following the EFAR catalogue, so extracted object N has the total errors of object N of the database, taking 
into account the correlated terms of the total errors (see Paper II) properly. Therefore, the error on cr, Mg2 and Mgb' for the 
extracted object N is computed considering the spectroscopic runs where object N of the database was observed, generating 
proper run correction errors and determining the final error taking into account the run weights. 

Finally, three additional sets of simulations (Cases 45, 46 and 47 of Table |l| are generated using uniform instead of 
gaussian distributions, and zero input peculiar velocity fields. Following the procedure described above, we consider the values 
of log R e , log a, (SB e ), (Ti, o"2, CT3 and of the coefficients of the FP as Case of Table |l|, and Eq. to specify the directions of 
the principal components of the distribution in the (log7? e , log a, (SB e )) space. For Case 45 we add gaussian errors generated 
as above, but we assume a uniform distribution of the FP parameters, extending ±v3cri, ±V3cr2, ±\/3cr3 in the wi,{?2,«3 
directions respectively from the mean values. This choice preserves the covariance matrix of the gaussian case. For Case 46 
we use gaussian distributions of the FP parameters, but measurements errors drawn using a uniform distribution extending 
±v3 the estimated rms. For Case 47 we use uniform distributions for both the FP parameters and the errors, with the y/3 
scaling as above. 



3.2 1-dimensional Distributions 



As a first example we discuss the use of the 1 Dimensional algorithm of §2.4. We generated 99 samples of the EFAR database 
(see Paper V) and considered the cr distribution. Fig. |l| shows the results as a function of the cut applied to the data. For cuts 
well below the true mean of the distribution, the ML algorithm with or without normalization correction estimates the mean 
and rms of the distribution with small bias. The simple mean, which does not take into account the selection weights of the 
data, slightly overestimates the true value, because objects with small a tend to have small selection probabilities. The simple 
rms, corrected for measurement errors as suggested by Arkitas and Bershadi (1996), underestimates slightly the true values, 
not taking into account the selection weights. As the cut is increased to values similar to the true mean of the distribution, 
the ML algorithm with normalization correction estimates the mean and rms with small bias, with the statistical errors (for 
the single sample \/99 larger than the error bars shown) becoming increasingly large. With cuts one sigma or more larger 
than the mean, the ML algorithm starts to fail, biasing the mean low and the sigma high. However, the bias is always smaller 
than the random errors. Both the ML algorithm without normalization correction and the simple estimates bias the mean to 
larger values and the rms to smaller values. 
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Figure 1. The results of the simulations of the <r distribution. The panels to the left show the mean (top panel) and rms (bottom panel), 
averaged over 99 simulations, as derived by the ML algorithm with normalization correction for a cut- The errorbars are l/%/99 smaller 
than the random error expected for one simulation. The central panels show the mean (top) and rms (bottom) as derived by the ML 
algorithm without normalization correction for cr cut . The panels to the right show the simple mean (top) and rms (bottom, corrected for 
measurement errors). The horizontal dotted lines give the input values. The vertical dotted lines mark the positions of the input mean 
value, and ±cr, ±2cr away from it. 



As it is well known (Press et al. 1986), the likelihood analysis offers in principle a simple recipe to compute confidence 
intervals for the fitted parameters. The na one- dimensional confidence interval is fixed by the contours of the likelihood 
function where: 

lnC max -\nC = 0.5n 2 . (29) 

Figure ^ shows four contours of constant likelihood in the (/i, a) plane for one of the simulations of Figure [l] analyzed by 
the ML algorithm with normalization correction a cu t = 100,150,200,250 km/s respectively. The thick errorbars mark the 
position of n and a, averaged over the 99 simulations of Fig. El, and show their rms. A straight application of Eq. H would 
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produce n — 1 one- dimensional confidence intervals for fx and a a factor ss 2 smaller than the rms measured by the simulations. 
Inspection of Eq. reveals that since the mean weighting of each galaxy is (w) = EitUi/JV, the na one- dimensional confidence 
interval is reached at a 0.5(w)n 2 distance from the maximum likelihood value. This factor is taken into account in Fig. [jj 
where the scaled n = 1 one-dimensional confidence intervals match approximately the rms measured by the simulations. The 
two-dimensional la confidence region, set by A In C = 1.15 (to) , is slightly larger. However, Fig. [l3| shows that the mean weight 
fluctuates in the range 1.5 — 2, making the estimate of confidence intervals from the likelihood function uncertain for a given 
dataset. In addition, this analysis does not allow to study the small residual biases in the ML estimates. Therefore, in the 
following we will estimate confidence intervals and residual biases using Monte Carlo simulations only. 



3.3 The Mg 2 -Mgb' relation 

The simplest test of the 2-dimensional algorithm is performed considering the Mg2-Mg&' relation. We generated 99 samples 
without taking into account the EFAR selection function. The input parameters are the ones derived in Paper II. Two cases 
are examined, one with a realistic error distribution and one with a factor 5 smaller errors. Figure |^(d) shows one of the 
realizations with realistic errors. Figure |(g) shows the results averaged over the 99 simulations. In both the cases of small 
and realistic errors the ML algorithm retrieves the input parameters with small biases and with small errors. The Y — X, 
X — Y and bisector regressions fail to derive the correct slopes and zeropoints, because of the finite orthogonal scatter and 
the non- negligible errors. When the method of Akritas and Bershady (1996) is used to correct for the errors (dotted cross), 
the Y — X result is less biased, but with large random errors. Note that the method cannot be used when selection weights 



are present. The orthogonal residual fit gives a slightly biased answer as expected from Eq. B6 , with a nearly perfect result 
when the errors are small. As expected from Eq. |l2|, the errors on a and b are perfectly anticorrelated. 

Figs. [|(a), (b), (c), (e) and (f) show the histograms of the results of the ML parameters (for the case of realistic errors). 
Not only are the slope and the mean values determined with high precision, but so also is the intrinsic scatter (orthogonal 
and parallel to the relation), despite the non-negligible measurement errors (the observed scatter orthogonal to the relation 
is more than a factor two larger than the intrinsic one). 



3.4 The Mg 2 -cr and Mg6' -a relations 

The tests on the Mg2-cr and Mg&'-cr relations were performed using the EFAR mock catalogues with realistic error distri- 
butions. The input parameters are the ones derived in Paper V. Figure ^(a) shows one of these Mg2-cr realizations. Data at 
small a or Mg2 have larger selection weights. Fig. ^(b) shows that the spread in the errors is large. The distribution of the 
selection probabilities (Fig. §(c)) has a peak at the completeness value (S = 1) with a prominent tail down to the selection 
probability cut (S = 0.1). The effective number of datapoints per selection weig ht bin (Fig. §(d)) is approximately constant. 
The simulations (full lines) reproduce closely the observed dataset (dotted lines) analyzed in Paper V. 

Fig. H(e) shows the slopes and zero-points derived from the various regressions and the ML algorithm, considering points 



with selection probabilities larger than 0.1 for the Mg2-cr relation. As already seen in 3. J, the regressions give biased results 
due to the errors (as discussed in § ^.l| ) , while the ML algorithm is nearly free of biases and accurate. 

Figure |^ illustrates the effects of the selection weighting and of the cut in a. When the analysis is performed including 
all the available datapoints to the largest selection weights and without any likelihood cut (open circle), the parallel and 
orthogonal scatter are derived best. However, the rms on all the parameters are rather large, because the presence of some 
datapoints with very high weights can bias the analysis of some realizations. It is therefore prudent to limit the analysis to 
datapoints with weights not larger than 10 (e.g., Si > 0.1) and clip points with low likelihood (— InC > 0, triangular stars). 
This allows us to determine the slope and the zero-point of the relation with small biases and with small variance. The price 
is a slight underestimate of the orthogonal and parallel spread, <ji and <72. Fitting points with weights not larger than 5 (filled 
triangles) and the same likelihood clipping, one observes that the mean quantities Mg 2 and log a are biased to systematically 
larger values, because the small galaxies have on average larger selection weights (see Fig. ^). The largest bias in the mean 
quantities is obtained weighting all datapoints equally, independently of the selection (open square). These tests justify the 
procedure adopted in Paper V to study the Mg2-cr and Mg6'-<r relations, where we consider galaxies with Si > 0.1 and apply 
likelihood clipping. 

To illustrate the effect of the cut in a, we examine the extreme case where only galaxies with a > 180 km/s are considered 
in the analysis. When no correction is applied (filled square), the slope and the zeropoint of the relation are slightly biased, 
the mean quantities log a and Mg 2 are biased to larger values, while the parallel spread is biased low (similar to what is seen 
in Figure [I]). The biases are fully corrected for when we apply the normalization corrections appropriate for a C ut = 180 km/s 
(open star), with a slight increase of the variance. The cut applied to the EFAR database (a > 100 km/s) is in any case low 
compared to the mean value of the sample. Finally, plot (d) of Figure |^ shows the mean values and rms (over the simulations) 
of the mean (over the datapoints) logarithmic likelihood and selection probability. The cross shows the values derived for the 
EFAR sample examined in Paper V, with errorbars equal to their rms divided by the square root of number of the galaxies 
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Figure 2. The contours of constant likelihood in the (fj,, a) plane for one of the simulations of Figure |l| analyzed by the ML algorithm 
with normalization correction a au t = 100, 150, 200, 250 km/s respectively. The projections (shown as dotted lines for the a C ut = 100 case) 
of the four full-line countours give the one-dimensional 1, 2, 3 and 4 sigma confidence intervals. The dotted contours show the 1 sigma 
confidence regions in two dimensions. The small crosses mark the maximum likelihood solution. The applied mean selection weights (w) 
scaling factors are given in each panel (see text). The dashed lines show the input values of the mean /i and rms a. The dots show the 
positions of maximum likelihood solutions of the simulations of Figure |l]. The thick errorbars mark the position of fi and cr, averaged 
over the 99 simulations of Figure |l], and show their rms. 

considered. The simulations match the mean likelihood of the EFAR sample within the estimated errors, justifying a posteriori 
the assumption of a gaussian distribution. The mean selection probability of the simulations (see Figure ^) is only slightly 
larger than the ones of the real sample. 

Similar results are obtained when considering the Mgb'-a relation. Figure ^(f) shows again that the linear regressions give 
biased results. The Y — X and orthogonal fitting are very similar, because of the shallow slope. The ML with probability cut 
Si > 0.1 and likelihood clipping gives nearly unbiased results. Histograms similar to Fig. ^(a)-(f) for the Mg2-o" and Mgb'-a 
relations have been already presented in Paper V. 
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3.5 The Fundamental Plane 

In this section we describe the tests of the ML algorithm used to derive the FP solution of Paper VII and the peculiar velocities 
of the EFAR sample. In §3.5.1 we assess the superiority of the ML algorithm with respect to the linear regression methods, 
when the a and b parameters of the FP are determined. In §3.5.2 we study the precision of the FP parameters derived using 
the ML algorithm, justifying the stn 
distribution of galaxies in the (log R e 
adopted in P 
data sample. 



the ML algorithm, justifying the strategy adopted in Paper VII 

log a, (SB C )) space 

adopted in Paper VII. Finally, 



3.5.3 justifies a posteriori the assumption of a gaussian 



3.5.4 tests the bias correction scheme for the peculiar velocities 



3.5.5 investigates how well coherent and random motions can be measured from the EFAR 



3.5.1 Method comparison 

As we have done above, we start comparing the performances of the "classical approaches" (linear regressions and orthogonal 
fitting) and the ML algorithm, focusing on mock catalogues generated using parameters close to those derived for the EFAR 
data sample (see Paper VII and Case of Table 0). We consider the 29 best clusters of the EFAR sample (see Paper VII) 
and we do not add peculiar velocities (i.e., the input peculiar velocity field is zero, Si npu t = 0). In each of the four methods 
considered we solve simultaneously for the coefficients a and b (see Eq. ^) and the vector Sj. The shifts 5j are determined 
by assuming that the mean of Sj over the 29 clusters is zero, which fixes the value of the parameter c. Objects with a < 100 
km/s, or with selection diameter Dw less than 12.6 kpc (i.e., logD^cut = 1.1), or with selection probabilities less than 0.1, 
are excluded from the fits. Therefore, this procedure mimics closely the "fiducial" solution of Paper VII ("Case 1" of Table 
4). 

Fig. ^ shows the mean values and rms (over 99 simulations) of the Fundamental Plane coefficients a and b as determined 
using regressions on logi? e , logo - , (SB e ), orthogonal fitting and ML. As can be expected from the previous discussion, linear 
regressions and orthogonal fitting are inadequate to determine the coefficients of the Fundamental Plane. All five methods 
determine b quite accurately, with the (SB e ) and logi? e regressions giving the largest systematic deviations (ss ±0.02). Only 
the ML algorithm, however, is able to estimate a with a bias smaller than the expected statistical uncertainty. As expected 
from Appendix H|, the logi? e and {SB e ) regressions and (less severely) the orthogonal fitting underestimate a, while the 
logo - regression overestimates it. As pointed out by Isobe et al. (1990), the logi? e regression gives the smallest rms. The ML 
algorithm determines a to 0.09 and b to 0.013 statistical accuracy. 

Figure shows how well the five methods measure the input peculiar velocity shifts Sj (all equal to zero) of the EFAR 



clusters. Similarly to the procedure described in §3.5.4 for the ML algorithm, the parameters a, b and c of the FP determined 
from the regressions on the subset of 29 best clusters for each separate simulation are used to compute the peculiar velocities 
of the remaining clusters. 

Let us first focus on the clusters with log D^ j — Sw,j < 1.18. These are the clusters that statistically have at least one 
galaxy with log Dw ~ 1.1 and selection probability larger than 0.1. In this case, the biasing influence of errors and selection 
effects is small for all the methods. The absolute values of the mean Sj are typically smaller than 0.02 dex. The rms over 
the EFAR clusters of the mean Sj over the 99 simulations is w 0.008. The log R e regression gives slightly poorer results. The 
statistical errors on the single cluster measurements are in any case larger, in the range 0.02-0.04 dex for the best-populated 
clusters, and up to 0.15 dex for the poorer populated clusters. The mean over the 99 simulations of the standard deviation 
(over the EFAR clusters) of the Sj values differs for the various methods. The log R e regression gives the smallest value (0.043 
dex, or a typical 10.4% accuracy on the determination of the distance of a single clusters), similar to the (SB e ) regression 
(0.045 dex, 10.9% distance error). The ML algorithm and the orthogonal regression perform only slightly worse (0.049 dex, 
or 11.9% distance error). The logcr regression gives the largest (0.065 dex, or 16% distance error) random error. 

The ML algorithm (and in order of increasingly importance the orthogonal, (SB e ) and log7? e regressions) derives Sj 
biased to low values when logD^ — 8w,j > 1.18. The clusters where this happens are so distant that their smallest galaxies 
in the EFAR sample have log Dw > 1-1. Fig. ^ shows that these are clusters where we sample less than half of the galaxy 



distribution. From what discussed in §3.2 we expect the ML algorithm to progressively fail in this regime. In §3.5.4 we correct 
for this residual bias due to the difference in the way the FP galaxy distribution is sampled in different clusters through 
simulations. In any case, the standard deviations (over the EFAR clusters) of the Sj values do not change, remaining w 4 
times larger that the typical systematic errors. 

Therefore we conclude that overall the ML method provides the best solution, keeping low the systematic and random 
errors on both the FP parameters and the Sj values. In addition, it is more robust against outliers (since we can identify and 
remove the extreme low-likelihood objects, see Paper VII). 
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3.5.2 Variant cases and systematic errors 

In this section we study the precision of the FP parameters derived using the ML algorithm, to justify the strategy adopted 
in Paper VII. Table |f| summarizes the tests performed, giving the mean values (over 99 simulations) of the number of clusters 
and galaxies (rounded to integer numbers) and the parameters of the Fundamental Plane derived for 28 cases. Table |^ lists 
the corresponding rms. Figure [] shows the results graphically. 

Case gives the parameters of the Fundamental Plane as derived in Paper VII for the standard fit. These are the input 
parameters used to generate the simulations considered here. Case 1 is the standard fit discussed above. Fig. 2 of Paper VII 
presents the histograms of 1000 simulations for this case. The mean values derived from the 99 simulations used here do not 
differ by more than 1.5% from those derived from the 1000 simulations. As discussed in Paper VII, there are small residual 
biases in the fitted parameters: a is biased low by 6% (i.e., since the mean value of a recovered from the simulations is smaller 
than the input value, we infer that the true value of a for the EFAR sample must be larger by ~ 6% than what found in 
Case 0), b is biased low by 2%; c is biased high by 4%; logT?,., log a and (SB e ) are all biased high, by 0.036 dex, 0.07 dex 
and 0.05 mag respectively; the scatter about the FP <ti is under-estimated by 0.006 dex, or 1.4%; the widths 02 and 0-3 of the 
galaxy distribution in the FP are biased by 0.049 (low) and 0.009 (high) dex respectively. These biases are all less than or 
comparable to the rms width of the distribution, so that although they are statistically significant (i.e. much greater than the 
standard error in the mean), they do not dominate the random error in the fitted parameters. Therefore in Paper VII we do 
not correct for these biases, since they are small and have negligible impact on the derived distances and peculiar velocities 



(see §3.5.4) 



Case 2 is similar to Case 1: we do not add peculiar velocities and we do not fit objects with a < 100 km/s, or with 
selection diameter Dw less than 12.6 kpc (i.e., logDwmt = 1-1), or with selection probabilities less than 0.1. However, in each 
simulation the list of fitted clusters is not fixed to the 29 clusters considered in Case 1, rather it is restricted to the clusters 
that have 6 or more galaxies. The number of fitted clusters fluctuates from 20 to 29 from simulation to simulation. The mean 
FP parameters do not show differences from Case 1. 

Cases 3 to 8 examine fits performed on the subset of the clusters having 10 or more galaxies (Case 3) to clusters having 
3 or more galaxies (Case 8). The mean FP parameters derived for Case 3 are the ones with the smallest biases. However, 
only 5-6 clusters for ~ 70 — 80 galaxies are fit. As a consequence, the statistical error on the parameters is nearly a factor two 
larger than for Case 1. On the contrary, Case 8 fits a larger number of galaxies and clusters than Case 1, and therefore gives 
slightly smaller (~ 20%) statistical errors. However, the mean values of the FP parameters are biased more than for Case 1. 
In particular, the value of o\ is spuriously small, as offsetting the FP with spurious peculiar velocity suppresses the apparent 
scatter. We conclude that the list of clusters adopted in Paper VII is a reasonable compromise between the need of obtaining 
bias-free parameters and minimizing the statistical errors. 

Cases 9 to 12 explore the effects of different logDwcut- Low logDvvcut give slightly smaller a coefficients; in addition, 
log R e is too large, since essentially no normalization correction is applied. However, the statistical errors are smaller, because 
the sample size is maximized. In contrast, large logDwcut give larger a coefficients and too small log_R e . In addition, the 
statistical errors are larger, because the sample size is reduced. We conclude that again Case 1 is the reasonable compromise 
both in terms of bias and statistical errors. This choice, however, forces us to apply additional corrections when deriving the 



peculiar velocities of the whole cluster sample (see ^3.5.4) 



Case 13 ignores selection probabilities altogether and applies a uniform weight to all galaxies, resulting in an effective 
over- weighting of the larger galaxies. As a consequence, the mean values of log7i e , log a, and (SB e ) are biased to higher 
values. A similar effect, but of reduced amplitude, is observed for Case 14, where galaxies with selection probabilities lower 
than 0.2 are excluded. On the contrary, Case 15 includes essentially all galaxies in the fit, irrespective of their selection 
probabilities. The biases in the recovered parameters are reduced, but the statistical errors increase, because of the presence 
of a few (deviant) points with large weights. To conclude, weighting the data points with their inverse selection probability, 
and excluding galaxies with selection probabilities lower than 0.1 is a reasonable compromise that minimizes the biases and 
the statistical errors of the recovered parameters. 

Case 16 shows that excluding the (few) galaxies with large errors on a does not affect the results of the fits, both the 
mean and rms values of the parameters. 

Case 17 excludes not only the galaxies rejected from the standard fit, but also galaxies with low likelihoods (In C < 0); this 
results in a highly biased fit, with low a coefficient, an artificially lowered FP scatter and a substantially narrower distribution 
in the FP. Case 18 replaces the individual error estimates for all measured quantities with uniform (average) errors. This has 
little effect on the derived parameters, but the statistical errors are larger. 

Case 19 allows an extra degree of freedom by permitting the orientation of the major axis of the galaxy distribution 
within the FP to be fitted, rather than specified a priori. The vectors V2 and C3 of Eq. |27| read now: 

VI = fl + U2,2^2 + "2,3^3, (30) 
#3 = "3 + V3 2X2 + X3. 
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The only parameters affected are 02 and 03 which are biased high and low respectively. Fig. |l(j shows the histogram of 
the recovered second and third component of vector V2, and first and second component of vector V3. The input parameters 
(derived from the input a and b values and Eq. ^) are recovered with residual biases that are smaller than the statistical 
error. Therefore in Paper VII we use the simplifying approximation of Eq. ^ 

Cases 20 to 23 explore the effects of errors on the parameters of the selection functions determined in Paper I. Case 20 
considers mock catalogues where the selection probabilities used in the fitting procedure have been computed adding +0.1 
dex to the true log Dw- As a consequence, the selection weights are systematically higher than the true ones (see Fig. |l3| ) 
and the mean number of selected galaxies (i.e. with selection probabilities larger than 0.1) is systematically smaller (see Table 
In Case 21 we added -0.1 dex to the true log Dw, getting selection weights systematically smaller than the true ones (see 
Fig. |l^) and therefore a larger mean number of selected galaxies. Cases 22 and 23 use selection weights distorted by adding 
+0.1 dex (Case 22) and -0.1 dex (Case 23) to the selection probability widths 8w The largest systematic effects are observed 
for Case 20, where systematically smaller values of a, b, log a, and systematically larger values of c and log R e (because more 
small galaxies are excluded from the fit) are obtained. See §3.5.3 and Fig. |l^ for a discussion of the effects on the distribution 
of likelihood probabilities. 

Cases 24 to 28 refer to mock catalogues of the EFAR sample with peculiar velocity fields. Case 24 has a random peculiar 
velocity field (always the same for the performed 99 simulations) of 0.05 dex rms amplitude. We considered this particular 
field (see Fig. [H], where we plot the peculiar velocities of all clusters, the 29 FP plus the remaining ones) because its mean 
bulk motions of the clusters at positive (the Hercules Corona Borealis sample, HCB) and negative Galactic longitudes (the 
Pisces Perseus Cetus sample, PPC) are slightly non-zero and of opposite sign, reminiscent of what measured in Paper VII. In 
addition, its mean velocities projected along the LP and SMAC dipole directions (see Fig. [l^ for the peculiar velocity field 
of all clusters) are also similar to what seen in Paper VII. Case 25 has a random field (always the same for the performed 99 
simulations) of the same rms amplitude as Case 24 (0.05 dex) plus the bulk flow motions determined by Lauer and Postman 
(1994, hereafter LP, see Fig. [l^ for the peculiar velocity field of all clusters). LP measured the bulk motion of the Abell Clusters 
with cz < 15000 km/s with respect to the Cosmic Microwave Background to be 764 km/s in the direction of I = 341, b = 44 
(reanalysis by Colless 1995), using the brightest cluster galaxies L — a relation. The resulting peculiar velocity field has a 
prominent asymmetry of the mean bulk motions of the HCB and PPC clusters (see Fig. |l2] ), which gives an apparently larger 
dipole motion in the LP direction (see Fig. |l^). Case 26 adds on top of a random field of 0.05 rms amplitude (always the same 
for the performed 99 simulations) the bulk motion determined by Hudson et al. (1999, hereafter SMAC). SMAC determined 
the motions of the Abell Clusters with cz < 12000 km/s using the FP, finding Vbulk = 630 km/s towards / = 260, b — — 1. 
The resulting peculiar velocity field has small mean bulk motions of the HCB and PPC clusters (see Fig. Il3 for the peculiar 
velocity field of all clusters), and an approximately zero dipole motion in the SMAC direction (see Fig. [18|). Case 27 has a 
pure LP bulk motion without a random component. Case 28 has a pure SMAC bulk motion without a random component. 
In all five cases the FP parameters are recovered without additional biases with respect to Case 1. The statistical errors are 
slightly larger, especially for the log R e parameter. The mean number of galaxies varies slightly, as galaxies near the selection 
cut might get included or not in the sample accordingly to the positive or negative peculiar velocity of the clusters they belong 
to (see Sect. [j.l| ). We conclude that the ML algorithm recovers the FP parameters correctly also when peculiar velocity fields 
are present. 

The last three sets of simulations (Cases 45, 46 and 47) refer to input distributions uniform in the FP parameters, in the 
errors, or both (see Sect, [[j]), with zero peculiar velocity fields. Inspection of Tables [l] and ^| and Fig. |{] shows that also in 
these cases the input parameters are recovered with small biases. The largest effects are seen for 02, where the ML gaussian 
algorithm overestimates the effects of the uniform error distribution, biasing 02 low. The rms values are similar to the ones 
derived for the previous cases: da, db, dc are marginally larger, doi marginally smaller. 

Summarizing, with only a few exceptions (noted above), the systematic differences in the fits derived for different cases 
are comparable to the random errors in the determination of the parameters for the standard case. Therefore, the uncertainties 
in our best-fit FP parameters are dominated by the random errors and not by systematic effects from the fitting method. 
In addition, the systematic effects explored in the simulations (Tables [l] and ^) reproduce, at least qualitatively, the trends 
observed in the fits to the actual EFAR data (Table 4 of Paper VII). 



3.5.3 The likelihood distribution 

We now investigate whether the assumption of a gaussian distribution of galaxies in the FP space is compatible with the 
EFAR data sample. Fig. [l^ (left panel) shows the cumulative distributions of the mean (i.e. normalized to the number of 
galaxies fitted) likelihoods of the simulations of Cases 1 to 8, 18 to 28, and 45 to 47 (top row, only galaxies of clusters used to 
determine the FP parameters) and Cases 1 and 20 to 47 (bottom row, all galaxies of the clusters for which p eculiar velocities 
have been computed). Note that Cases 29 to 44, not listed in Tables ^ and ^, are discussed in detail in § 3.5.4 They are 
perturbed versions of Case 1. 

Except for Cases 20 and 21, all simulations (including the uniform distribution cases 45, 46 and 47) give similar cumulative 
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Table 1. The mean values (over 99 simulations) of the parameters of the Fundamental Plane derived for various cases. 



Case 


N cl 


N S al 




a 




6 


c 


log/? e 


logo- 


(SB e ) 








02 




C3 


Notes 





29 


255 


1 


.223 


0. 


336 


-8.666 





.770 


2 


.304 


19 


.71 





.064 


1 


.995 





.610 


standard fit Paper VII 


1 


29 


240 


1. 


.150 


0. 


329 


-8.360 


0. 


.808 


2 


.311 


19 


.77 





.058 


1 


.946 





.643 


standard fit 1000 simulations 


1 


29 


240 


1. 


.168 


0. 


330 


-8.418 


0. 


.799 


2 


.310 


19 


.76 





.058 


1 


.950 
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distributions of the mean likelihoods. Depending on the exact case, a mean likelihood equal to or larger than that of the EFAR 
sample is observed in up to 30% of the simulations. Case 45 (a uniform distribution of the FP parameters and gaussian errors) 
gives only slightly smaller mean likelihoods. Case 20 gives systematically higher mean likelihoods, and Case 21 systematically 
lower, than that of the EFAR sample. Inspection of the central and left panel plots of Fig. [l^ clarifies that the cumulative 
distribution of the mean unweighted likelihoods are similar for all the simulations with gaussian error distributions, predicting 
that values larger than that of the EFAR sample are observed in 70 to 95% of the cases. The Cases 46 and 47, where the error 
distributions are uniform, produce systematically 20-30% larger values, because the ML gaussian algorithm is overstressing 
the role of errors. As noticed before (see Fig. |E]), the mean selection weights of all simulations, except Case 20, are smaller 
than that of the EFAR sample. 

Therefore, we conclude that the gaussian modeling is a reasonable description of the distribution of the EFAR galaxies 
in the (logi? e , logo - , (SB e )) space and of the error distribution. The EFAR data set does not allow to discriminate between 
a gaussian or a uniform distribution of the FP parameters. Error distributions with tails slightly stronger than gaussians 
are hinted. Values of log D^y slightly larger (by ~ 0.04 dex) than that measured in Paper I are needed to match the mean 
selection weight of the EFAR sample. 



3.5.4 Peculiar velocities and bias corrections 

Once the parameters of the best-fitting gaussian distribution are determined using a subset of clusters (hereafter, the FP 
clusters), the ML algorithm is used with these parameters fixed, to compute the peculiar velocities of the remaining clusters. 
Tables ^[ ^ and ^ give the average residuals (over 99 simulations) from the input peculiar velocities < AS > =< S — Si npu t > 
for a number of cases. Table [] considers Cases 1 and 20 to 28 of Table [l], listing also the values of < AS > given in Table 7 of 
Paper VII, based on 1000 simulations of Case 1, and the statistical precision reached for Case 1 and 99 simulations. Tables ^ 
and |E| give the results for 19 additional Cases, from 29 to 47. The 16 Cases from 29 to 44 are obtained by perturbing one of 
the FP parameters by plus or minus one sigma, keeping the others fixed to the best solution value. Cases 29 and 30 perturb 



© 0000 RAS, MNRAS 000, 000-000 



Peculiar Motions of Early-Type Galaxies. VI. 15 



Table 2. The rms values (over 99 simulations) of the parameters of the Fundamental Plane derived for various cases. 



Case dN c i dN ga i da db dc dlog/2 e dlogcr d(SB e ) do\ d<J2 do"3 Notes 

1 4 0.089 0.013 0.354 0.025 0.016 0.060 0.006 0.158 0.067 standard fit 1000 simulations 

1 4 0.089 0.013 0.347 0.027 0.018 0.048 0.006 0.168 0.066 standard fit 99 simulations 

2 1 10 0.093 0.014 0.360 0.027 0.018 0.048 0.006 0.172 0.068 uses clusters with N gat > 6 

3 1 8 0.141 0.021 0.533 0.047 0.028 0.087 0.011 0.246 0.108 uses clusters with N gat > 10 

4 1 10 0.117 0.016 0.445 0.038 0.025 0.063 0.008 0.198 0.089 uses clusters with N gal > 8 

5 1 9 0.106 0.015 0.411 0.033 0.022 0.057 0.007 0.179 0.079 uses clusters with N gal > 7 

6 2 9 0.085 0.013 0.338 0.026 0.017 0.044 0.006 0.156 0.062 uses clusters with N gal > 5 

7 1 7 0.079 0.012 0.328 0.022 0.014 0.043 0.005 0.143 0.054 uses clusters with N gal > 4 

8 2 6 0.074 0.012 0.309 0.022 0.014 0.041 0.005 0.136 0.054 uses clusters with N gai > 3 

9 3 0.067 0.011 0.274 0.019 0.009 0.048 0.005 0.162 0.039 no D w cut is applied 

10 3 0.067 0.011 0.275 0.018 0.009 0.046 0.005 0.162 0.039 D Wcut = 6.3 kpc 

11 5 0.097 0.014 0.380 0.041 0.028 0.048 0.007 0.167 0.075 D Wcut = 14.1 kpc 

12 6 0.097 0.014 0.366 0.068 0.055 0.059 0.007 0.157 0.109 D Wcut = 15.9 kpc 

13 4 0.068 0.011 0.281 0.016 0.010 0.037 0.005 0.122 0.048 uses no selection weighting 

14 5 0.077 0.013 0.323 0.021 0.012 0.046 0.006 0.149 0.057 uses galaxies with Si > 0.2 

15 3 0.101 0.014 0.385 0.047 0.032 0.064 0.006 0.179 0.104 uses galaxies with Si > 0.001 

16 3 0.093 0.013 0.362 0.029 0.019 0.048 0.006 0.171 0.070 uses galaxies with 8a < 0.1 

17 10 0.150 0.026 0.667 0.030 0.017 0.064 0.010 0.268 0.080 excludes galaxies with In C< 

18 4 0.106 0.014 0.385 0.031 0.024 0.047 0.009 0.175 0.114 uses uniform errors for all galaxies 

19 4 0.087 0.013 0.333 0.025 0.017 0.049 0.006 0.161 0.040 also fit third axis of FP 

20 5 0.082 0.014 0.359 0.022 0.014 0.054 0.006 0.172 0.059 all log D° shifted by +0.1 

21 3 0.086 0.012 0.325 0.027 0.018 0.042 0.006 0.145 0.069 all log D w shifted by -0.1 

22 3 0.087 0.012 0.329 0.027 0.018 0.041 0.006 0.146 0.067 all S w shifted by +0.1 

23 5 0.086 0.015 0.379 0.026 0.016 0.056 0.006 0.174 0.064 all 5 W shifted by -0.1 

24 4 0.087 0.014 0.357 0.032 0.020 0.059 0.006 0.127 0.068 random peculiar velocity field 

25 4 0.096 0.014 0.410 0.030 0.019 0.058 0.006 0.156 0.072 LP dipole plus random component 

26 4 0.078 0.013 0.354 0.029 0.016 0.063 0.005 0.139 0.065 SMAC dipole plus random component 

27 3 0.087 0.016 0.398 0.028 0.019 0.057 0.006 0.170 0.070 pure LP dipole 

28 3 0.096 0.012 0.378 0.029 0.019 0.056 0.005 0.147 0.079 pure SMAC dipole 

45 4 0.096 0.011 0.356 0.027 0.020 0.045 0.005 0.126 0.072 FP uniform distribution 

46 3 0.093 0.016 0.402 0.022 0.015 0.045 0.006 0.138 0.058 Error uniform distribution 

47 4 0.100 0.017 0.438 0.031 0.021 0.053 0.005 0.116 0.073 FP and Error uniform distribution 



the a coefficient, 31 and 32 b, 33 and 34 logi? e , 35 and 36 logo - , 37 and 38 (SB e ), 39 and 40 en, 41 and 42 02, 43 and 44 03, 
respectively. Cases 45, 46 and 47 refer to the simulations with uniform distributions in the FP parameters, the errors or both. 

Fig. ^ shows < AS >-< < AS > > as a function of logD w . Here < < AS > > is the average residual over the FP 
clusters. It is zero by construction for Cases 1, 20 to 28 and 45 to 47. The general trend of increasingly negative residuals at 
large values of logD w is easily explained. The clusters at higher redshifts (large logD w ) lack small galaxies and therefore 
have mean values of log R e larger than log R e , or 8 — log R e — log i? e progressively more negative. 

Once the zeropoint offset (given by the residual < AS >averaged over the FP clusters) is subtracted off, the simulations 
29 to 44 give residuals compatible with Case 1 within the statistical errors. Therefore, the residual bias corrections are 
robust against small errors in the FP solution. The largest variations of the zeropoint are observed when the values of 
log _R e , log a and (SB e ) are perturbed. As expected, one finds < < AS > >~ ±0.025 = ±dlogi? e , when logi? e is perturbed, 
< < AS > >Ki ±0.02 = +a dlog a, when log a is perturbed, and < < AS > ±0.02 = +6 d{SB e ) when (SB e ) is perturbed. 
Errors in log D w affect the selection weighting scheme and therefore the determination of < AS > for the clusters (at higher 
redshifts) that have the largest log D w . In some cases the errors in < AS > are larger than the statistical precision of Case 
1. Positive errors (+D W , Case 20) produce too negative < AS >, since more small galaxies than with the correct values of 
log D^y are excluded from the fits. The contrary happens for negative errors (— Dyy, Case 21). Errors in Sw affect also < AS >, 
especially at higher redshifts. In this case negative errors (-dD w , Case 23) produce too negative < AS >, because a sharper 
(i.e., a smaller Sw) selection function eliminates more small galaxies. Positive errors (dD w , Case 22) produce the reverse. 

The presence of a peculiar velocity field (Cases 24 to 28) can also affect < AS >, because redshift distances are used to 
convert the measured selection cuts from arcsec to kpc. In particular, clusters sampled by a small number of galaxies are prone 
to small number statistics fluctuations, as galaxies near the selection cutoff might not be always included in the simulation 
catalogues. 

Uniform distributions of the FP parameters, of the errors or both do not affect < AS > systematically, within the 
statistical errors. 

Finally, Fig. [| illustrates the effect of the log Dwcut value on the required peculiar velocities. When no log Dwcut is 
applied (Case 9), the residual biases at high logD^ are larger than when logDwcut is equal to 1.1 (Case 1). The biases are 
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reduced if a larger logDwcut (1-2, Case 12) is chosen, but the statistical errors increase. In addition, as discussed in §3.5.2, 
this choice gives larger systematic and statistical errors for the FP parameters. 

We conclude that the peculiar velocity biases are much smaller than the random errors in the peculiar velocities. The 
uncertainties in the peculiar velocity bias corrections as computed from the 1000 simulations of Case 1 used in Paper VII do 
not increase the random errors in the peculiar velocities. However, residual systematic effects at the level of 0.01 dex might 
still be present, especially for the clusters at higher redshifts, due to the uncertainties in the determination of the selection 
function. 



3. 5. 5 The estimation of the mean motions 

Given the large mean distance of the EFAR clusters, and the precision of the FP distance estimator (20% for the single galaxy, 
see Paper VII), the present dataset does not allow the determination of the peculiar motions of single clusters. Nevertheless, 
the mean motions of the EFAR clusters can be constrained. In the following we investigate this point using the subsample of 
clusters identified in Paper VII for the study of the bulk motions. To this belong the 50 EFAR clusters with 3 or more galaxies 
(hereafter Peculiar Velocity, PV, clusters), redshifts less than 15000 km/s, and peculiar velocity errors less than 1800 km/s. 



We compute the average mean velocities without applying the bias corrections investigated in § 3.5.4 (although this is done 
in Paper VII), to stress that our results are not affected much by these residual systematic uncertainties, and that random 
errors dominate the noise. 

We begin by examining the mean peculiar velocities < V >, < Vhcb > and < Vppc > of the sample. These are the 
mean over the PV clusters, over the PV clusters with positive (the Hercules Corona Borealis, HCB, sample), and over the 
ones with negative Galactic longitude (the Perseus Pisces Cetus, PPC, sample). Fig. |l6| shows that these average quantities 
are recovered with statistical errors of about 200 km/s. In general these errors are larger than the systematic differences from 
the input values. In particular, errors on the selection function parameters do not affect the result much (Cases 20 to 23) and 
uniform distributions of the FP parameters, of the errors or both do not change this conclusion (Cases 45 to 47). The largest 
systematic difference is observed for Cases 24 to 26 and Vppc, where the input values are ill-defined, since the rms of the 
cluster peculiar velocities is large (see Figs. [TT] and [t2| ). 

As a second test we consider the Supergalactic X, Y and Z components of the mean peculiar velocities in radial shells. 
Fig. [l7] shows the averages of 99 simulations obtained for Case 1 and 24 of Table [jj. Again, the input values are recovered with 
systematic differences smaller than the statistical errors. 

We determine now how well the EFAR cluster subsample can constrain the bulk flow motions determined by Lauer and 
Postman (1994) and Hudson et al. (1999). We fitted the peculiar velocities of the 50 clusters used in Paper VII with a bulk 
flow model Vbuik cos 8, where 6 is the angle between the dipole and cluster direction, determining Vbuik in a least squares sense. 
Fig. 15 of Paper VII shows these fits to the EFAR data. Fig. ^ shows that when simulations with no peculiar velocities are 
considered (Cases 1 and 20 to 23 of Table [j]), the recovered mean Vbuik is compatible with zero, with statistical errors of about 
160 km/s. Errors on the selection function parameters do not affect the result much (Cases 20 to 23), as uniform distributions 
of the FP parameters, of the errors or both (Cases 45 to 47). Similarly, when the simulations with pure LP or SMAC dipole 
motions are considered (Cases 27 and 28 respectively) , the recovered mean Vbuik is consistent with the input values within the 
statistical uncertainties. Fig. 16 of Paper VII shows the histograms of the recovered bulk flow amplitudes for 500 simulations. 
Only one out of these simulations of the LP flow, and none of SMAC yields a Vbuik less than the value measured from the 
EFAR data in Paper VII. In contrast, the simulations with a random peculiar velocity component (Case 24) give mean Vbuik 
values similar to the observed values, both for the LP and the SMAC directions. Simulations with a random component on 
top of the LP and SMAC bulk flows (Cases 25 and 26) produce large mean Vbuik along the LP direction, and little net motion 
along the SMAC direction. 

Finally, we investigate whether not only the mean motions, but also their global intrinsic fluctuations could be determined 



with the EFAR cluster sample. To this purpose we use the one-dimensional ML algorithm of §2.4: we consider the peculiar 
velocities derived for the 99 simulations of Case 24 of Table ^ and their projections along the supergalactic X, Y and Z 
coordinates with their estimated errors, separately for the whole sample of PV clusters, and the northern HCB and southern 
PPC subsets. For each of these datasets the ML algorithm estimates the intrinsic (gaussian) dispersion. Fig. [H)] shows the 
results. The derived rms are systematically underestimated, up to more than a factor two, when the input values are of the 
order of ~ 600 km/s, and by ~ 25% if larger input values are considered. We conclude that cluster sample is too small to 
allow an accurate determination of these quantities. 

To summarize, the EFAR cluster sample allows to measure the mean motions of the HCB and PPC regions with 200 
km/s precision, and to constrain the bulk flow along the LP dipole direction strongly. It is less suited to constrain the SMAC 
result and does not allow the accurate determination of the random peculiar motion component. 
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Figure 3. The results of the Mg2-Mgb' simulations. The panels show the histograms of (a) the slopes, the mean values of (b) < Mgb 1 > 
and (c) < Mg2 >, the orthogonal and parallel spread (e) a\ and (f) <T2 as derived from the ML algorithm. Input values, means and 
rms of the 99 simulations are also given. Panel (d) shows one of the 99 Mg2-Mgb' realizations. Panel (g) shows the values of the slopes a 
and zeropoints b, averaged over 99 simulations and derived using the Y-X (y), the X-Y (x), bisector (b) and orthogonal (o) regressions, 
plus the ML algorithm (m). The errorbars show the random errors expected for one simulation. The dotted cross shows the values and 
the random errors of the slopes a and zeropoints b, averaged over 99 simulations, derived using the Y-X regression, taking into account 
errors as in Akritas and Berschady (1996). The horizontal and vertical dashed lines show the input values. The other dotted lines of plot 
(g) terminating to the small crosses show the values obtained when the errors are reduced by a factor 5. The crosses for the ML and the 
orthogonal regression case overlap. 
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Figure 4. The results of the Mg2-cr and Mgfe' —c simulations. Panel (a) shows one of the 99 Mg2— ^realizations. The size of the points 
is proportional to their selection weight 1/S. Panel (b) shows the distribution of the Mg2 errors. Panel (c) plots the distribution of 
selection probabilities S for the simulation shown in (a) and (b) (full line) and for the EFAR database (dotted line) analyzed in Paper 
V. Panel (d) plots the effective number of data points N/ S as a function of S for the simulation shown in (a) and (b) (full line) and for 
the EFAR database (dotted line) analyzed in Paper V. Panel (e) and (f) show the values of the slopes a and zeropoints (b), averaged 
over 99 simulations and derived using the Y-X (y), the X-Y (x), bisector (b) and orthogonal (o) regressions, plus the ML algorithm (m), 
for the Mg2-cr and Mgfc'— a relations, respectively. The errorbars show the random errors expected for one simulation. The horizontal 
and vertical dotted lines show the input values. 
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Figure 5. The effects of the selection weighting and of the cut in a. The plots show the mean values (over 99 simulations) of (a) the 
slope and zeropoint, (b) the mean values, (c) the orthogonal and parallel scatter as determined by the ML algorithm, (d) the mean 
logarithmic likelihood and mean selection probability. The errorbars show the random errors expected for one simulation. The different 
types of points refer to the cases discussed in the text. The cross in panel (d) shows the mean values derived for the EFAR sample 
examined in Paper V, with errorbars equal to their rms, divided by the square root of number of the galaxies considered. 
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Figure 6. The parameters of the FP relation. The plot shows the mean coefficients < a > and < b > (over 99 simulations) derived 
from the linear regressions on logi? e , logcr, (SB e ), the orthogonal minimization and the ML algorithm. The errorbars show the rms. The 
dotted lines show the input values. 
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< 1.18. The input peculiar velocity shifts are all zero. 
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Figure 8. One of the simulations of §3.5.1 Panel (a) shows the correlation between log R e -0.3{SB e ) and logDvy. The full line shows 
the input mean value log R e -0.3{SB e ) , the dashed dashed line the value derived by the ML algorithm. The dotted lines correspond to 
the given values of log D\y cu t . Panel (b) shows the values of the selection parameters log ^ and Sw,j • The 29 best clusters are shown 
circled. Panel (c) shows the cumulative distribution of log R e -0.3{SB e ) values. Panel (d) shows the cumulative distribution of logD^-. 
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Figure 9. The fitted FP parameters for each case in Table |l| showing the distributions and correlations for various pairs of parameters. 
Each case is numbered as in the Table. The dotted lines show the input parameters. The errorbars show the rms over 1000 simulations 
derived for Case 1. 
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Fig ure 10. The distributions of the second and third component of vector V2, and first and second component of vector 113 (see Eq. 
pot) resulting from fitting 99 simulations allowing the orientation of the major axis of the galaxy distribution within the FP to be a free 
parameter (Case 19). The input parameters of the simulations are given at the head of each panel (and indicated by the vertical dotted 
line), followed by the mean and rms of the fits to the simulations. 
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Figure 11. The input peculiar velocity field of the simulations of Case 24 of Table [l| Each cluster is identified by its Cluster Assignment 
Number (hereafter CAN, see Paper II). Clusters in the southern Galactic hemisphere are plotted at negative rcdshifts. The resulting 
mean velocities (with statistical errors) for the whole sample, the clusters at positive (the Hercules Corona Borealis sample, HCB) and 
negative Galactic longitudes (the Pisces Perseus Cetus sample, PPC) are given with their rms. 
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Figure 12. The input peculiar velocity fields of the simulations of Cases 25 to 28 of Table [lj Labels as in Fig. [llj 
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Figure 13. The cumulative distributions of the mean likelihoods (left), mean unweighted likelihoods (center), mean selection weights 
(right). The top row refers to the galaxies of the clusters used in the determination of the FP parameters, the bottom row to all galaxies 
of the clusters for which peculiar velocities are computed. The vertical dashed lines show the values of EFAR sample, as determined in 
Paper VII. The numbers identify the Cases discussed in the text. 
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Figure 14. The residuals selection bias determined from 99 simulations and Cases 1 and 20 to 28 of Tables [j] and|i| and Cases 29 to 44 
of Tables ^ and |i| The top panel shows the residuals < A<5 >=< S — &i npu t > from the input peculiar velocity field 5i nput (8i nput = 
for Cases 1, 20 to 23, and 29 to 44; see text, Table [l] and Figs. |i"l||i2| for Cases 23 to 28), with the average residual over the FP clusters 
subtracted off, as a function of logD^,. Each Case is identified by its number (see Tables [l| [j| and ^). The error bars show the statistical 
precision reached for Case 1 with 99 simulations. The bottom panel shows the average residual measured for the FP clusters as a function 
of the Case number. The labels provide an additional mnemonic identifier. 
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Figure 15. The residual selection bias < AS > — << AS » determined for Case 9 (no logDivcut applied, top) and Case 12 
(log Dnr cut =l.2, bottom), as a function of log D^. The error bars show the statistical precision reached with 99 simulations. The crosses 
show the mean values for Case 1. 
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Figure 16. The simulations of the mean peculiar velocities < V > (mean velocity of the sample of 50 Clusters defined in Paper VII, 
top panel), < Vhcb > (mean velocity of the clusters with positive Galactic longitude, the HCB sample, middle panel) and < Vppc > 
(mean velocity of the clusters with negative Galactic longitude, the PPC sample, bottom panel). The points show the mean values (over 
99 simulations) and their rms for Cases 1, 20 to 28, and 45 to 47 of Table M. The dotted bars show the input values. 



© 0000 RAS, MNRAS 000, 000-000 



Peculiar Motions of Early-Type Galaxies. VI. 31 



_ All 



500 



A 
> 
V 



-500 - £ 



1 1 r 

■ SGX 
• SGY 

A SGZ 



HCB 



5 3 12 9 7 9 5 



~i 1 1 r 



1 1 



If 



7 6 4 4 



— i 1 1 1 1 r 

PPC 



■u - 



4 2 5 3 3 5 3 



J_ 



_i i i_ 



10000 15000 10000 15000 10000 15000 

cz (km/s) cz (km/s) cz (km/s) 



_ All 
1000 - 



* o - 



A 
> 

V 



-1000 



+ 1 



SGX 

SGY 
-SGZ 



HCB 



PPC 



_l I L 



T 



10000 15000 10000 15000 10000 15000 

cz (km/s) cz (km/s) cz (km/s) 



Figure 17. The simulations of the mean peculiar velocities in radial shells. The clusters are grouped in 7 redshift ranges: the first is 
4000-8000 km/s, the next five cover 8000 km/s to 13000 km/s in 1000 km/s steps, and the last is 13000-15000 km/s. The left panel 
shows the whole sample of 50 PV clusters, the middle panel shows the HCB clusters, and the right panel shows the PPC clusters. The 
points show the mean values (over 99 simulations) of the peculiar velocities and their rms for Cases 1 (top) and 24 (bottom) of Tabic [l| 
The Supergalactic X, Y and Z components are shown as filled squares, circles and triangles respectively (with small offsets in redshift 
for clarity). In the bottom panels the dotted, short dashed and long dashed lines show the input values. The number of clusters in each 
redshift range is indicated at the bottom of each top panel. 
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Figure 18. The simulations of the bulk flow motions of Lauer and Postman (1994, top) and Hudson et al. (1998, bottom). The points 
show the mean values (over 99 simulations) of Vtulk derived using the 50 clusters used in Paper VII for the determination of the EFAR 
mean peculiar velocities, for 11 Cases identified by the numbers in Table |l| The errorbars show the rms. The dotted lines show the values 
of Vf,ulk as derived by LP and SMAC. The full lines show the values of Vf, u ;fc as derived in Paper VII. 
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Figure 19. The histograms of the peculiar velocity dispersions derived from the one-dimensional ML analysis of the 99 simulations of 
Case 24 of Table [l| The dispersions are shown separately for the HCB, PPC and PV cluster samples, for the radial peculiar velocities 
and their projections along the supcrgalactic X, Y and Z projections. The labels give the input dispersions, the results averaged over the 
99 simulations and their rms. 
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Table 3. The mean residual selection bias < A5 > for Cases 1 (for 1000 and 99 simulations, with the rms error of the 99 simulations), 
and 20 to 28 of Table [j] 
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0. 
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0. 
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-0 
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-0 


.0053 





.0016 


-0 
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0. 
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0. 
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.0212 


-0 
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.0016 


29 
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0. 


.0066 
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.0081 


-0 


.0073 


-0. 


.0069 


-0 
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-0 


.0129 


-0 
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-0 


.0149 


-0 


.0042 


-0 


.0017 


30 


-0. 


.0135 


-0. 
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-0 


.0241 


-0 


.0140 


-0. 


.0158 


-0 


.0214 
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-0 
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.0197 


-0 
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-0 


.0219 


-0. 
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-0 
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-0 
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-0 
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.0059 


0. 


0024 





.0022 


-0 


.0029 
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-0. 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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0. 
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0. 
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-0 
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-0 
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-0 
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.0147 
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-0 
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-0 


.0247 


-0 


.0517 


-0 


.0342 


-0 


.0271 


-0 


.0349 


45 


+0. 
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-0 
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-0 
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-0 


.0329 


0. 


.0072 


-0. 


.0449 


-0 


.0327 
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.0338 


-0 


.0410 


-0 
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-0 


.0502 


-0 
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.0346 


-0 
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.0011 
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-0 


.0165 
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.0034 


-0. 


.0173 


-0. 


.0053 


-0. 


.0063 


-0 


.0134 


-0 


.0129 


-0 


.0190 


-0 


.0168 


-0 


.0109 


-0 


.0071 


59 


-0 


.0047 


-0. 


.0021 





.0031 


-0 


.0011 


-0. 


.0025 


-0. 


.0039 





.0018 





.0013 


-0 


.0107 


-0 


.0089 


-0 


.0072 


-0 


.0012 


60 


-0 


.0251 


-0. 


.0155 


0. 


.0039 


-0. 


.0233 


-0 


.0126 


-0. 


.0139 


-0 


.0197 


-0 


.0179 


-0 


.0243 


-0 


.0210 


-0 


.0165 


-0 


.0268 


61 


+0. 


.0083 


0. 


0062 


0. 


.0045 


0. 


.0099 


0. 


0084 


0. 


.0074 





.0088 





.0013 


-0 


.0016 





.0025 


-0 


.0041 


-0 
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-0 


.0066 
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0027 





.0044 


-0 


.0015 


-0 
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-0 
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-0 


.0053 


63 


+0. 


.0010 


-0. 


.0077 





.0105 


-0. 


.0055 


-0. 


.0042 


-0. 


.0065 


-0 


.0042 


-0 


.0192 


-0 


.0096 


-0 


.0192 


-0 


.0044 


-0 


.0146 


64 


-0 


.0309 


-0. 


.0218 


0. 


.0064 


-0. 


.0266 


-0 


.0150 


-0 


.0182 


-0 


.0221 


-0 


.0304 


-0 


.0425 


-0 


.0439 


-0 


.0306 


-0 


.0383 


65 


-0 


.0015 


0. 


0080 


0. 


.0028 


0. 


.0130 


0. 


0085 


0. 
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.0061 


-0 
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-0 
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-0 
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Table 3. Continued. 



CAN 


1 (1000) 


1 (991 




d<5 




20 




21 




22 




23 




24 




25 




26 




27 




28 
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0. 
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.0070 
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.0006 


-0 
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-0 
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0040 
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-0 
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-0 
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-0 
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-0 
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-0 
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68 
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-0 
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-0 
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-0 
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-0 
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-0 
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69 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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.0356 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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-0 
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+0.0026 
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0. 


0156 
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.0038 





.0042 


-0 
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-0 
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.0042 


-0 


.0141 


-0. 
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-0 
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-0 


.0111 


-0 
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79 


-0.0022 


0.0017 


0. 


0037 





0071 





0014 





.0022 





.0030 


-0 


.0030 


-0 


.0037 


-0 


.0147 


-0 


.0012 


-0 


.0032 
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-0 
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.0022 





.0087 


-0 
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.0012 





.0058 


-0 
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-0 
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-0 


.0047 
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83 


-0.0024 


0.0012 





.0029 


0. 
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0. 


0012 





.0016 





.0039 


-0 


.0013 


-0 


.0049 


-0 


.0151 





.0003 


-0 


.0004 


90 


+0.0035 


0.0071 





0020 





0117 





0074 





.0068 





.0115 





.0037 


-0 


.0002 


-0 


.0062 





.0062 





.0038 
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Table 4. The mean residual selection bias < A8 > for Cases 29 - 38. 



VN 




29 




30 




31 




32 




33 






+da 




-da 




+dft 




-db 


+dlog R e 


1 





.0059 


-0. 


.0027 


0. 


.0037 


-0. 


0003 


0. 


.0267 


2 


-0 


.0320 


-0 


.0422 


-0 


.0341 


-0 


.0401 


-0 


.0121 


3 


-0 


.0108 


-0 


0225 


-0. 


.0149 


-0 


0183 


0. 


.0084 


4 


-0 


.0052 


-0 


.0139 


-0 


.0077 


-0 


0113 


0. 


0154 


5 


-0 


.0091 


-0. 


.0209 


-0 


0119 


-0 


.0184 


0. 


.0108 


6 


0. 


.0062 


0. 


0008 


0. 


.0047 


0. 


.0024 


0. 


.0285 


7 


-0 


.0118 


-0 


.0232 


-0 


.0151 


-0 


.0197 





.0076 


8 





0065 


0. 


0002 


0. 


0059 


0. 


.0010 


0. 


.0284 


9 


0. 


.0095 





.0075 


0. 


.0077 


0. 


.0094 


0. 


0343 


10 





.0043 


-0. 


.0008 


0. 


.0031 


0. 


0006 


0. 


.0268 


11 


-0 


.0044 


-0 


.0102 


-0 


.0054 


-0 


.0092 


0. 


.0178 


12 





.0063 





.0010 





.0045 





.0029 


0. 


.0287 


13 





.0091 





.0030 


0. 


.0069 





.0053 


0. 


.0311 


14 


-0 


.0002 


-0 


0029 


-0 


0008 


-0. 


0022 


0. 


.0235 


15 





0137 


0. 


.0086 


0. 


0127 


0. 


.0097 


0. 


.0362 


16 


0. 


.0048 





0009 





.0032 


0. 


.0027 


0. 


.0279 


17 


0. 


0011 


-0. 


0056 


-0 


.0010 


-0 


0033 


0. 


.0228 


18 


-0 


.0012 


-0. 


.0053 


-0. 


.0028 


-0. 


.0036 


0. 


0218 


19 


0. 


.0019 


-0 


.0009 


0. 


.0018 


-0 


.0007 


0. 


.0274 


20 


0. 


.0089 





.0053 





.0077 





.0066 





.0321 


21 


0. 


.0057 


0. 


.0029 





.0049 


0. 


.0037 


0. 


.0293 


22 


0. 


.0049 


0. 


0022 





.0047 





.0025 


0. 


.0286 


23 


0. 


.0047 


-0. 


.0008 


0. 


.0030 


0. 


.0009 


0. 


.0269 


24 





.0042 


-0 


.0012 


0. 


.0025 





.0006 


0. 


.0265 


25 


-0 


.0054 


-0 


.0124 


-0 


.0072 


-0 


0106 


0. 


.0161 


26 


-0 


.0035 


-0. 


.0109 


-0. 


.0061 


-0 


.0083 


0. 


0178 


27 





.0149 





0119 


0. 


.0130 


0. 


.0137 


0. 


.0383 


29 


-0 


.0045 


-0 


.0121 


-0 


.0050 


-0 


.0114 





.0168 


30 


-0 


.0133 


-0. 


.0215 


-0 


.0155 


-0 


.0193 


0. 


.0084 


31 





.0109 


0. 


.0082 


0. 


.0097 


0. 


.0094 


0. 


.0348 


32 





.0048 


-0 


0013 





.0030 





.0006 


0. 


.0268 


33 


-0 


.0021 


-0 


.0076 


-0. 


.0032 


-0 


.0066 


0. 


0201 


34 


-0 


.0045 


-0 


0136 


-0 


.0066 


-0 


.0114 


0. 


0160 


35 


0. 


.0053 


-0. 


.0006 


0. 


.0034 


0. 


0015 


0. 


.0274 


36 


-0 


.0045 


-0 


0123 


-0. 


.0067 


-0. 


0101 


0. 


0166 


37 





.0070 





.0027 





.0060 





.0038 


0. 


.0299 


38 





.0064 


-0 


.0003 


0. 


.0046 





.0015 


0. 


.0281 


39 





.0125 


0. 


.0078 


0. 


0117 


0. 


0088 


0. 


.0352 


40 





.0022 


-0 


.0026 


-0 


.0003 


-0 


.0001 





.0248 


42 


-0 


.0018 


-0 


.0082 


-0 


.0036 


-0 


.0063 


0. 


.0200 


43 


-0 


.0151 


-0. 


.0264 


-0 


0178 


-0 


.0237 


0. 


.0043 


44 


-0 


.0274 


-0. 


.0409 


-0 


.0307 


-0 


.0376 


-0. 


.0091 


45 


0. 


.0131 





0080 





.0118 





.0093 


0. 


.0356 


46 





.0018 


-0 


.0046 


-0 


.0004 


-0 


.0023 





.0236 


48 


0. 


.0035 


-0 


.0030 





.0017 


-0 


.0011 


0. 


.0253 


49 


-0 


.0009 


-0 


0094 


-0 


.0027 


-0 


.0074 


0. 


.0199 


50 





.0096 


0. 


.0045 


0. 


.0079 


0. 


.0062 


0. 


0321 


51 


0. 


.0109 


0. 


0044 


0. 


.0085 





.0067 


0. 


.0329 


52 


0. 


.0192 





.0135 





0171 


0. 


.0155 


0. 


0413 


53 


-0 


.0173 


-0 


0290 


-0 


.0209 


-0 


.0254 


0. 


.0019 


55 





.0026 





.0001 


0. 


.0032 


-0 


0005 


0. 


.0264 


56 


-0 


.0280 


-0 


.0378 


-0 


.0310 


-0 


.0348 


-0. 


.0079 


57 


0. 


.0048 


-0. 


0015 


0. 


.0030 


0. 


.0003 


0. 


.0267 


58 


-0 


.0037 


-0. 


0115 


-0. 


.0053 


-0 


.0098 


0. 


0175 


59 





0010 


-0. 


.0054 


-0. 


.0008 


-0. 


.0034 


0. 


.0229 


60 


-0 


.0107 


-0 


.0204 


-0. 


.0124 


-0. 


.0185 


0. 


.0095 


61 





.0075 


0. 


.0049 





.0064 





.0060 


0. 


.0312 


62 





.0075 


0. 


.0004 


0. 


.0054 


0. 


.0026 


0. 


.0290 


63 


-0 


.0060 


-0 


.0093 


-0. 


.0061 


-0. 


.0092 


0. 


.0186 


64 


-0 


.0173 


-0. 


.0264 


-0. 


0194 


-0. 


.0242 


0. 


.0032 


65 





.0109 





0050 


0. 


.0088 





.0072 


0. 


.0330 





34 




35 




36 




37 




38 


dlog R e 


+d(SB e ) 


-d(SBe) 


+dlog a 


-dlog a 


-0 


.0233 


-0 


0176 





.0209 


-0 


,0103 





,0125 


-0 


.0621 


-0. 


.0564 


-0 


.0177 


-0 


,0510 


-0 


.0225 


-0 


.0416 


-0. 


.0360 





.0028 


-0 


,0327 


-0 


.0001 


-0 


.0345 


-0 


.0288 





.0097 


-0 


,0221 





.0026 


-0 


.0413 


-0 


.0354 





.0049 


-0 


.0199 


-0 


,0104 


-0 


.0214 


-0. 


.0159 





.0230 


-0 


.0134 





.0206 


-0 


.0424 


-0. 


.0367 





.0020 


-0 


.0328 


-0 


.0014 


-0 


.0216 


-0 


.0161 





.0229 


-0 


.0108 





.0190 


-0 


.0173 


-0 


.0114 





.0284 


-0 


.0078 





.0228 


-0 


.0232 


-0 


.0176 





.0212 


-0 


.0152 





,0197 


-0 


.0322 


-0. 


.0268 





.0123 


-0 


.0223 





.0076 


-0 


.0213 


-0. 


.0157 





,0232 


-0 


,0134 





.0212 


-0 


.0189 


-0. 


.0132 





.0255 


-0 


.0079 





.0208 


-0 


.0265 


-0 


.0208 





.0177 


-0 


.0172 





.0146 


-0 


.0138 


-0 


.0081 





.0305 


-0 


.0016 





.0233 


-0 


.0221 


-0. 


.0164 





.0223 


-0 


.0127 





.0188 


-0 


.0272 


-0 


.0214 





.0171 


-0 


.0159 





.0110 


-0 


.0282 


-0 


.0226 





.0162 


-0 


.0192 





,0133 


-0 


.0264 


-0. 


.0203 





.0213 


-0 


.0159 





,0191 


-0 


.0179 


-0. 


.0122 





,0264 


-0 


.0098 





.0245 


-0 


.0207 


-0 


.0151 





.0237 


-0 


.0119 





.0204 


-0 


.0214 


-0 


0157 





.0229 


-0 


,0128 





,0191 


-0 


.0231 


-0 


0175 





.0214 


-0 


,0151 





,0187 


-0 


.0235 


-0. 


.0179 





.0210 


-0 


,0155 





,0181 


-0 


0339 


-0. 


.0283 





.0105 


-0 


,0241 





.0064 


-0 


.0322 


-0 


.0266 





.0122 


-0 


,0221 





.0078 


-0 


.0117 


-0. 


.0060 





.0326 





.0005 





.0259 


-0 


0333 


-0. 


.0275 





.0110 


-0 


.0215 





.0063 


-0 


.0432 


-0 


.0373 





.0025 


-0 


.0338 





.0001 


-0 


.0157 


-0 


.0103 





.0294 


-0 


.0075 





,0275 


-0 


.0232 


-0. 


.0176 





.0212 


-0 


.0153 





.0190 


-0 


.0299 


-0. 


.0243 





.0146 


-0 


.0200 





.0088 


-0 


.0340 


-0. 


.0283 





.0102 


-0 


,0231 





.0036 


-0 


.0226 


-0 


.0170 





.0219 


-0 


,0143 





.0189 


-0 


.0334 


-0 


.0278 





.0110 


-0 


,0247 





.0073 


-0 


.0201 


-0. 


.0145 





.0243 


-0 


,0116 





.0213 


-0 


.0219 


-0. 


.0162 





.0224 


-0 


,0121 





.0184 


-0 


.0148 


-0 


.0091 





.0296 


-0 


.0043 





,0252 


-0 


.0252 


-0. 


.0194 





.0191 


-0 


.0163 





,0149 


-0 


.0300 


-0. 


.0245 





.0146 


-0 


.0202 





,0101 


-0 


.0457 


-0 


.0399 


-0 


.0015 


-0 


.0332 


-0 


,0085 


-0 


.0591 


-0 


0533 


-0 


.0148 


-0 


.0482 


-0 


,0190 


-0 


.0144 


-0. 


.0089 





.0301 


-0 


.0057 





.0262 


-0 


.0263 


-0. 


.0206 





.0179 


-0 


.0150 





.0130 


-0 


.0247 


-0. 


0190 





.0195 


-0 


.0134 





.0150 


-0 


.0301 


-0 


.0244 





.0142 


-0 


.0205 





,0119 


-0 


.0179 


-0 


.0123 





.0264 


-0 


.0085 





,0226 


-0 


.0176 


-0. 


.0120 





.0272 


-0 


.0076 





,0222 


-0 


.0087 


-0 


.0032 





.0359 





.0028 





,0293 


-0 


.0481 


-0 


.0425 


-0 


.0037 


-0 


.0368 


-0 


,0092 


-0 


.0236 


-0. 


.0183 





.0210 


-0 


.0155 





.0195 


-0 


.0579 


-0. 


.0521 


-0 


.0136 


-0 


.0452 


-0 


.0202 


-0 


.0234 


-0 


.0176 





,0210 


-0 


.0121 





,0148 


-0 


.0325 


-0 


.0269 





.0118 


-0 


.0226 





,0079 


-0 


.0271 


-0 


.0214 





.0172 


-0 


,0179 





.0136 


-0 


.0405 


-0. 


.0348 





.0038 


-0 


.0306 





.0006 


-0 


.0188 


-0. 


.0133 





,0257 


-0 


.0106 





.0215 


-0 


.0210 


-0. 


.0153 





,0233 


-0 


.0107 





.0180 


-0 


.0340 


-0 


0281 





,0126 


-0 


.0237 





.0094 


-0 


.0468 


-0 


0411 


-0 


,0026 


-0 


.0358 


-0 


.0080 


-0 


.0170 


-0. 


0113 





,0273 


-0 


.0073 





.0241 
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Table 4. Continued. 







29 




30 




31 


32 


33 


34 


35 


36 


37 


38 






+da 




-da 




+db 


-db 












dlogcr 


66 





0095 


0. 


.0054 


0. 


0082 


0.0068 


0.0325 


-0.0175 


-0.0120 


0.0270 


-0.0086 


0.0238 


67 


-0 


.0027 


-0 


0112 


-0. 


.0055 


-0.0083 


0.0181 


-0.0319 


-0.0263 


0.0125 


-0.0223 


0.0081 


68 





.0014 


-0 


.0031 


-0 


.0000 


-0.0016 


0.0241 


-0.0258 


-0.0202 


0.0185 


-0.0140 


0.0119 


69 


-0 


.0152 


-0. 


0310 


-0. 


0204 


-0.0257 


0.0054 


-0.0515 


-0.0454 


-0.0007 


-0.0400 


-0.0060 


70 





0042 


-0. 


0002 


0. 


0033 


0.0009 


0.0271 


-0.0229 


-0.0173 


0.0214 


-0.0130 


0.0168 


71 





.0170 





.0121 





.0157 


0.0135 


0.0399 


-0.0106 


-0.0051 


0.0344 


0.0001 


0.0291 


72 


-0 


.0048 


-0 


.0123 


-0. 


.0070 


-0.0101 


0.0164 


-0.0336 


-0.0280 


0.0109 


-0.0206 


0.0054 


73 





.0131 


0. 


.0078 





.0106 


0.0104 


0.0357 


-0.0148 


-0.0090 


0.0299 


-0.0043 


0.0257 


74 


-0 


.0310 


-0. 


0398 


-0. 


0334 


-0.0374 


-0.0102 


-0.0607 


-0.0549 


-0.0160 


-0.0494 


-0.0226 


75 


-0 


.0017 


-0 


.0068 


-0. 


.0034 


-0.0051 


0.0208 


-0.0292 


-0.0235 


0.0150 


-0.0174 


0.0093 


76 


-0 


.0216 


-0. 


.0308 


-0 


.0237 


-0.0287 


0.0004 


-0.0528 


-0.0467 


-0.0057 


-0.0420 


-0.0102 


77 


0. 


0160 


0. 


.0116 


0. 


0138 


0.0138 


0.0388 


-0.0112 


-0.0056 


0.0333 


-0.0034 


0.0323 


78 


-0 


.0053 


-0. 


0105 


-0. 


.0073 


-0.0085 


0.0171 


-0.0329 


-0.0272 


0.0115 


-0.0241 


0.0084 


79 





.0044 


-0 


.0010 





.0020 


0.0014 


0.0267 


-0.0233 


-0.0177 


0.0211 


-0.0153 


0.0191 


80 





.0042 





0001 


0. 


.0029 


0.0016 


0.0272 


-0.0228 


-0.0172 


0.0216 


-0.0149 


0.0201 


82 





.0043 


-0 


.0000 


0. 


.0027 


0.0017 


0.0272 


-0.0228 


-0.0172 


0.0216 


-0.0144 


0.0190 


83 





.0032 


-0. 


.0009 


0. 


0019 


0.0006 


0.0262 


-0.0238 


-0.0181 


0.0205 


-0.0148 


0.0175 


90 


0. 


0087 


0. 


0054 


0. 


0074 


0.0069 


0.0321 


-0.0179 


-0.0122 


0.0264 


-0.0097 


0.0237 
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Table 5. The mean residual selection bias < 8 > for Cases 39-47. 



VN 




39 




40 




41 




42 




43 




44 




45 




46 




47 




+d<ri 




dcri 


+d<T2 




do- 2 


+dcr 3 




-d<T3 




UG 




GU 




uu 


1 


0. 


.0021 





0012 


0. 


.0018 


0. 


0015 





.0057 


-0 


.0031 


-0 


.0066 


-0 


.0131 


-0 


.0072 


2 


-0. 


.0377 


-0 


0365 


-0. 


.0368 


-0 


.0373 


-0 


.0334 


-0 


.0416 


-0 


.0250 


-0 


.0355 


-0 


.0236 


3 


-0 


.0172 


-0 


.0160 


-0. 


.0165 


-0 


.0167 


-0 


.0137 


-0 


.0203 


-0 


.0077 


-0 


.0115 


-0 


.0185 


4 


-0 


.0093 


-0 


0098 


-0. 


.0095 


-0 


0096 


-0 


.0059 


-0 


.0139 


-0 


.0071 


-0 


.0050 


-0 


.0022 


5 


-0 


.0146 


-0 


.0159 


-0. 


.0152 


-0 


.0154 


-0 


.0075 


-0 


.0227 





.0077 


-0 


.0026 


-0 


.0002 


6 





.0042 





0030 


0. 


.0036 


0. 


0035 





.0047 





.0020 





.0041 


-0 


.0027 





.0063 


7 


-0 


.0179 


-0 


.0169 


-0 


.0173 


-0 


.0175 


-0 


.0140 


-0 


.0217 


-0 


.0156 


-0 


.0138 


-0 


.0203 


8 


0. 


.0040 





.0029 





.0036 


0. 


.0032 





.0058 





.0004 


-0 


.0003 


-0 


.0037 


-0 


.0026 


9 


0. 


0103 





0069 


0. 


.0084 





.0086 





.0086 





.0083 





.0115 


-0 


.0003 





.0129 


10 


0. 


.0024 


0. 


.0013 


0. 


0019 





.0017 





.0029 





.0003 





.0042 





.0040 





.0068 


11 


-0 


.0066 


-0 


.0079 


-0 


.0071 


-0 


.0074 


-0 


.0054 


-0 


.0097 


-0 


.0111 


-0 


.0007 





.0035 


12 





.0044 





.0031 





.0037 





.0037 





.0048 





.0022 





.0069 





.0063 





.0080 


13 


0. 


.0070 





0053 


0. 


0061 





.0061 





.0087 





.0029 





.0045 





.0099 





.0108 


14 


-0 


.0003 


-0 


.0027 


-0. 


0015 


-0 


.0016 


-0 


.0008 


-0 


.0026 





.0086 





.0081 





.0052 


15 





.0123 


0. 


.0101 


0. 


.0113 





0110 





.0135 





.0083 





.0062 





.0074 





.0042 


16 


0. 


.0039 





.0021 





.0029 





.0029 





.0042 





.0013 





.0029 





.0031 





.0069 


17 


-0. 


.0016 


-0 


0027 


-0. 


.0022 


-0 


0022 





.0005 


-0 


.0055 


-0 


.0017 


-0 


.0065 





.0023 


18 


-0 


.0023 


-0 


0040 


-0. 


.0032 


-0 


0032 


-0 


.0021 


-0 


.0047 





.0038 





.0082 





.0110 


19 


0. 


.0019 


-0 


.0007 





.0006 


0. 


.0004 





.0011 


-0 


.0003 





.0232 





.0050 





.0141 


20 





.0081 





.0062 





.0071 





.0071 





.0080 





.0058 





.0037 





.0052 





.0046 


21 


0. 


.0056 


0. 


0032 


0. 


.0043 


0. 


0043 





.0050 





.0033 





.0126 





.0096 





.0070 


22 





.0047 





0025 


0. 


.0037 





.0035 





.0042 





.0027 





.0022 





.0066 





.0136 


23 





.0025 





.0015 





.0021 





.0018 





.0031 





.0003 





.0017 





.0025 


-0 


.0033 


24 


0. 


0021 





.0011 





0016 





.0015 





.0027 


-0 


.0001 


-0 


.0026 





.0035 





.0009 


25 


-0 


.0086 


-0 


.0091 


-0 


.0087 


-0 


.0090 


-0 


.0067 


-0 


.0117 


-0 


.0101 


-0 


.0099 


-0 


.0046 


26 


-0. 


.0068 


-0 


0076 


-0. 


.0072 


-0 


.0072 


-0 


.0051 


-0 


.0099 





.0012 





.0155 


-0 


.0006 


27 


0. 


0149 





0119 


0. 


0132 


0. 


.0134 





.0139 





.0125 





.0192 





.0217 





.0290 


29 


-0 


.0080 


-0 


.0085 


-0 


.0080 


-0 


.0086 


-0 


.0050 


-0 


.0121 





.0031 


-0 


.0047 


-0 


.0078 


30 


-0 


0173 


-0 


.0175 


-0 


.0173 


-0 


.0175 


-0 


.0151 


-0 


.0204 


-0 


.0145 


-0 


.0050 





.0020 


31 


0. 


0109 


0. 


0084 


0. 


.0095 


0. 


0096 





.0101 





.0088 





.0107 





.0030 





.0147 


32 





.0023 


0. 


.0014 


0. 


0019 


0. 


.0017 





.0030 





.0002 


-0 


.0029 





.0001 





.0058 


33 


-0 


.0045 


-0 


.0051 


-0 


.0048 


-0 


.0050 


-0 


.0032 


-0 


.0070 


-0 


.0051 





.0007 





.0016 


34 


-0 


0091 


-0 


.0089 


-0 


.0088 


-0 


.0092 


-0 


.0054 


-0 


.0134 


-0 


.0191 


-0 


.0027 


-0 


.0088 


35 


0. 


.0029 





0020 


0. 


.0025 


0. 


0024 





.0038 





.0006 


-0 


.0011 


-0 


.0023 





.0036 


36 


-0 


.0083 


-0 


0085 


-0. 


.0083 


-0 


0085 


-0 


.0066 


-0 


.0107 


-0 


.0077 


-0 


.0039 


-0 


.0033 


37 





.0058 





.0041 





.0049 





0048 





.0059 





.0036 





.0043 


-0 


.0018 





.0006 


38 





.0037 


0. 


.0025 


0. 


.0032 





.0029 





.0051 





.0004 





.0052 





.0066 


-0 


.0045 


39 


0. 


0112 


0. 


.0093 


0. 


0103 


0. 


0101 





.0119 





.0081 





.0139 





.0021 





.0042 


40 


0. 


.0006 


-0 


0009 


-0. 


.0003 


-0 


0001 





.0013 


-0 


.0022 





.0021 


-0 


.0023 





.0024 


42 


-0 


.0044 


-0 


.0055 


-0 


.0049 


-0 


.0050 


-0 


.0031 


-0 


.0073 





.0070 





.0053 


-0 


.0007 


43 


-0 


.0212 


-0 


.0202 


-0 


.0205 


-0 


.0210 


-0 


.0157 


-0 


.0266 


-0 


.0248 


-0 


.0258 


-0 


.0338 


44 


-0 


0353 


-0 


0329 


-0. 


.0338 


-0 


.0344 


-0 


.0290 


-0 


.0402 


-0 


.0308 


-0 


.0333 


-0 


.0279 


45 


0. 


0114 





0098 


0. 


0106 





.0105 





.0117 





.0090 





.0022 





.0018 





.0116 


46 


-0 


.0007 


-0 


.0020 


-0 


0013 


-0 


.0014 





.0012 


-0 


.0045 


-0 


.0037 





.0031 


-0 


.0038 


48 


0. 


.0012 


-0 


.0006 





.0003 


0. 


.0002 





.0029 


-0 


.0030 


-0 


.0013 


-0 


.0017 


-0 


.0052 


49 


-0 


0051 


-0 


.0051 


-0. 


.0049 


-0 


.0053 


-0 


.0025 


-0 


.0084 


-0 


.0077 


-0 


.0076 


-0 


.0032 


50 


0. 


.0079 


0. 


.0063 


0. 


.0071 





.0070 





.0087 





.0049 





.0078 





.0026 





.0032 


51 





.0083 





.0070 





.0077 


0. 


.0076 





.0094 





.0053 





.0087 





.0089 





.0002 


52 





.0175 





.0153 





.0163 





.0163 





.0181 





.0140 


-0 


.0004 





.0160 





.0217 


53 


-0 


.0240 


-0 


0223 


-0. 


.0230 


-0 


.0232 


-0 


.0187 


-0 


.0284 


-0 


.0216 


-0 


.0233 


-0 


.0286 


55 


0. 


.0025 


0. 


0004 


0. 


0014 





.0013 





.0020 





.0005 





.0045 





.0070 





.0057 


56 


-0 


0332 


-0 


.0325 


-0 


.0328 


-0 


.0329 


-0 


.0279 


-0 


.0387 


-0 


.0435 


-0 


.0214 


-0 


.0294 


57 


0. 


.0025 


0. 


.0009 





.0017 


0. 


0016 





.0042 


-0 


.0015 





.0029 





.0029 





.0087 


58 


-0. 


.0075 


-0. 


.0076 


-0. 


.0074 


-0. 


.0078 


-0 


.0050 


-0 


.0108 


-0 


.0091 


-0 


.0129 


-0 


.0035 


59 


-0. 


.0017 


-0. 


.0025 


-0. 


.0021 


-0. 


.0022 


-0 


.0002 


-0 


.0046 





.0012 


-0 


.0034 


-0 


.0018 


60 


-0 


.0158 


-0. 


.0153 


-0 


.0153 


-0. 


.0158 


-0 


.0124 


-0 


.0194 


-0 


.0189 


-0 


.0160 


-0 


.0153 


61 





.0074 





.0051 





.0062 





.0062 





.0068 





,0054 





,0186 





0001 





.0102 


62 


0. 


.0046 


0. 


.0035 


0. 


.0041 


0. 


.0039 





.0066 





.0008 


-0 


.0008 





.0053 





.0061 


63 


-0 


.0064 


-0. 


.0089 


-0. 


.0076 


-0. 


.0078 


-0 


.0071 


-0 


.0085 


-0 


,0132 





.0112 


-0 


.0030 


64 


-0 


.0218 


-0 


0219 


-0 


0216 


-0 


.0221 


-0 


.0185 


-0 


.0260 


-0 


.0264 


-0 


.0251 


-0 


.0169 


65 





0087 





.0074 





.0080 


0. 


.0079 





.0097 





.0058 





.0082 





.0020 





.0035 
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Table 5. Continued 



VN 




39 




40 




41 




42 


43 




44 




45 




46 




47 




+d<ri 




-dfjx 


+d(T2 




-d(T2 


^do"3 




-d(T3 




UG 




GU 




uu 


66 





.0084 


0. 


0067 


0. 


.0076 





0074 


0.0086 





.0060 





.0045 





.0064 





.0052 


67 


-0 


.0067 


-0. 


.0070 


-0. 


.0068 


-0 


.0070 


-0.0042 


-0 


.0103 


-0 


.0216 


-0 


.0057 


-0 


.0036 


68 


0. 


.0000 


-0 


.0016 


-0 


.0008 


-0 


.0009 


0.0008 


-0 


.0030 





.0025 





.0037 





.0033 


69 


-0 


.0239 


-0. 


.0223 


-0 


.0229 


-0 


.0231 


-0.0182 


-0 


.0291 





.0003 


-0 


.0312 


-0 


.0166 


70 





.0029 


0. 


0013 


0. 


0021 





0020 


0.0037 


-0 


.0000 


-0 


.0012 





.0113 





.0023 


71 


0. 


.0156 


0. 


.0137 


0. 


.0147 





.0145 


0.0161 





.0126 





.0194 





.0039 





.0136 


72 


-0 


.0080 


-0. 


.0091 


-0 


.0085 


-0 


.0086 


-0.0052 


-0 


.0126 


-0 


.0001 


-0 


.0072 


-0 


.0075 


73 





.0116 





.0095 





.0104 





.0105 


0.0120 





.0085 





.0099 





.0021 





.0008 


74 


-0. 


.0356 


-0. 


0352 


-0. 


0353 


-0 


0356 


-0.0321 


-0 


.0395 


-0 


.0398 


-0 


.0257 


-0 


.0391 


75 


-0 


.0035 


-0. 


.0049 


-0. 


.0042 


-0. 


.0042 


-0.0022 


-0 


.0067 





.0035 





.0057 





.0031 


76 


-0 


.0259 


-0 


.0265 


-0. 


.0260 


-0 


.0264 


-0.0238 


-0 


.0293 


-0 


.0233 


-0 


.0155 


-0 


.0208 


77 





.0148 


0. 


0129 


0. 


.0138 


0. 


0138 


0.0147 





.0126 





.0136 





.0057 





.0128 


78 


-0 


.0073 


-0. 


.0084 


-0 


.0078 


-0. 


.0079 


-0.0065 


-0 


.0097 


-0 


.0041 





.0037 





.0077 


79 


0. 


.0023 


0. 


0012 


0. 


0017 





0017 


0.0029 





0001 





.0079 





.0045 





.0029 


80 


0. 


.0030 


0. 


0015 





.0022 


0. 


.0022 


0.0031 





.0010 





.0054 





.0052 





.0055 


82 





.0030 





.0015 


0. 


.0022 





.0022 


0.0033 





.0006 





.0080 





.0039 





.0030 


83 





0021 





.0003 


0. 


0012 


0. 


0012 


0.0024 


-0 


.0005 





.0031 





.0082 





.0045 


90 





0083 


0. 


0060 





.0071 


0. 


.0071 


0.0079 





.0059 





.0069 





.0031 





.0088 
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4 CONCLUSIONS 

We have described the maximum likelihood gaussian algorithm we developed to investigate the correlations between the 
parameters of the EFAR database (the Mg2-Mg&' relation of Paper II, the Mgb'-cr and Mg2-<r relations of Paper V, the 
Fundamental Plane of Paper VII) and determine the cluster peculiar velocities of Paper VII. 

We performed extensive testing based on mock catalogues of the EFAR sample. We find that "canonical" methods based 
on a least-squares approach cannot cope with the challenge of a sample with a spread of a factor of two in redshift, with sizable 
selection effects, non-negligible and non-uniform errors, and explicit cuts. We quantify the size of the systematic biases these 
methods introduce. In contrast, the maximum likelihood gaussian algorithm takes into account errors, selection effects and 
the presence of explicit cuts, determining nearly-unbiased estimates of the slopes of the correlations and their intrinsic and 
parallel spread. Ten to thirty percent of the analyzed simulations have mean likelihoods larger than that of the EFAR sample, 
justifying the use of gaussian modeling. We derive the analytical solution of the maximum likelihood gaussian problem in N 
dimensions in the presence of small errors. We show that the residual systematic biases are always smaller than the statistical 
errors. We investigate in detail the effects of cluster sample selection, errors in the selection function parameters, and selection 
cuts. We explore the cases of uniform distributions of parameters and errors. We conclude that the mean peculiar motions of 
the EFAR clusters can be determined reliably. In particular, the large amplitude of the dipole motion measured by Lauer and 
Postman (1994) can be strongly constrained. 
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APPENDIX A: 1-DIMENSIONAL MODEL, LIMITING CASE 

Neglecting the presence of cuts, the maximization of Eq. ( ]l7| ) leads to: 

M= V 1 

2,2/ \2 
Oi + O- - (Xi - (J.) 



(Al) 
(A2) 



Here we solve Eqs. Al and A2 in the limiting case of small measurement errors [a%a). We get: 



a 



,2 

err 



fj, = X + (x - X er r)-^-, (A3) 



S 

V 2 = S 2 - olrr + 2cr 2 rr ^1 - ^f^J , (A4) 

where x = i^2i Xi /^ i )l (E V^O i s the selection weighted mean, x err = (E ■ Xi<7^/Si)/ a 2 /Si) is the error and selection 
weighted mean, a 2 rr = (E£r|/5i)/(Sl/Sj) is the selection weighted mean square error, s 2 = (EO 1 — x) 2 / <Si)/(Ei 1/Si) the 
selection weighted rms and s 2 rr = ( E o\{x — x) 2 / 'Si) / '(E» a 2 /Si) the selection and error weighted rms. Therefore the mean 
and rms are obtained by correcting the usual formulae with a term taking into account the spread in errors. 



APPENDIX B: THE 2-DIMENSIONAL MODEL, LIMITING CASE 

Here we present the solution which minimizes Eq. |l^ in the limiting case of small errors (with a diagonal error matrix), no 
peculiar velocities and no cuts, considering expansions to first order in the errors, for the two dimensional case. In Appendix 
|c| we derive the general result in the N dimensional case. To simplify the notation in the following we call the components 
(a;i,a;2) as (x,y)- For the mean values we find: 

x = *c°> + c&cv£V 0) - - <V 0) - yi^/A^, (Bi) 
y = y <0) + <x£,(-v£V 0) - 4 15 ) + vg\yW y«))/A< >, (B2) 

where i<°) = Q^.Xi/SJ/S, = {^yi/SJ/S, S = E, x£ ] = (E, z« = (£. x^J Si) / {So\ y ), 

<& = (Ei^i/Si)/^ <4 = (Ei°*/&ys> y* 1} = (EiW^i/St)/^), = (E^i/aVC^*). ^ = (E»(^- 

xWf/SMS, V 2 f = (E 4 (W - y (0) ) 2 /S0/S, 4 0) = (E^ - * (0) )fo - J/ (0) )/S0/S, and A (0) = ^»^» - V,f 2 . The 
covariance matrix V is: 

= VW-O*. (B3) 
+2aL(l - (V^al x>ex - V^a xy , ex )/A^), 

V 22 = vg+vv 

= V 2 f-a 2 y (B4) 
+24,(1 - (V^f'C - ^V IWV )/A(°)), 

^ + 

- a 2 ex (V 2 fa xy , ex ~V^a 2 yy . ex )/A^ (B5) 

- cr 2 (V (0) ct - V (0) cr 2 VA (0) 

= (Eifc* - ^ (0) ) 2 a 2 ,/S0/(5aL), <r 2 H , eH = (JJ.( W - yWyo&SiMSoSy), = - xC))^ - 

y )& x i/ Si)/{Sa 2 x ), <y xy , ey = (Ei^ ~~ * )(2/i ~~ y^) a li/Si)/(Sa 2 y ). As in the 1-dim case, we recognise the usual (er- 
ror corrected) zeroth-order plus a term taking into account the spread in the errors. With no spread in the errors we have 
V$ = and we recover the ansatz of Akritas and Bershady (1996). To first order Eq. ^|then reads: 

( T/ (i) T/ (i) t/(i) T/ (o) ,r(0) 

1+ ; V22 ~ Vl1 \^ V22 ~ Vl1 =\, IB(.) 
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where ai is Eq. evaluated using V^°K With no spread in the errors: 



(o) 



,(0)2 
'12 



Therefore underestimates the true slope if the error in the x direction is larger than the one in y. 



(B7) 



APPENDIX C: DERIVATION OF THE N-DIMENSIONAL CASE 

Let V be the intrinsic covariance matrix and A its inverse, Ei the error matrix of each point, Vi = V + Ei the observed 
covariance matrix of each point, with Ei « V. Then we have: 



Ai = (V + E,)' 1 = [Vil + V^Ei)]- 1 = (I + AE i )- 1 A~(I-AE i )A = A-AE i A, 
where the first order approximation (I + e) _1 ~ / — e was used. We get: 
detVi = detVdet(l + AE t ) = detV[l + tr(AEi)] = detV[l + E jk A jk Eij k ], 

where the first order approximation det(I + e) ~ f + tr(e) was used. Therefore to first order we get: 
detAi = detA[l ~ Ej k Aj k E i>jk ]- 
The equation for the likelihood reads: 
VdetA 



(2tt)"/ 2 



(1 - 0.5Ej k A jk E iJk ) exp 



-Xi (A — AEiA)xi/2 



l/Si 



(CI) 
(C2) 
(C3) 

(C4) 



where the n-dimensional m vectors of datapoints, and x the vector of the mean values. Taking the logarithm 

we get: 

1 



lnr = E£ 1 



St 



n I i 1 " t 

- ln(27r) + - IndetA - -T. jk A jk Ei, jk - -x t (A - AEiA)xi. 



The equation for the vector of the mean values reads: 
' " u L E,(A - AEiA){xi - S) k /S t = 0. 



dx k 

Setting x = s' ' + or 1 ', the zero-th order term is: 

s m = (E^ l /s I )/(sa/s I ). 

The first order term is then: 
- EiA^/Si = E ! A£ l A(f» - f (0) )/&, 
and therefore: 

= (E l £ s A(f (0) - xJ/SiJ/iZil/Si), 



which is equivalent to Eqs. B1-B2 for Ei diagonal and taking the zeroth order for A (see below). 
Let now consider the derivatives with respect to the A components: 



<91n£ 



<9A im 

Let first note that: 
f ddetA _ (-l) l+m detA (lim) 



1 ddetA 1^ Ajk la,T dA - , 1 - T . dAEjA > ) _ ,, 

2detA 8A lm 2^ jk * i < jk dA lm 2* 1 ( dA lm > Xl + 2** 1 0A !m jXl J ~ " ' 



A- 1 



det.4 <M im detA 

where ^4(;, m ) is the matrix obtained eliminating the row I and the column m. Since A is symmetric we get: 

detA dAim 

if I 7^ m, and: 

t ddetA _ 
detA dAu ~ 11 

otherwise. Similarly, we get 8Aj k /dAi m = 5jiS km + 5j m 8 k i if I ^ m and dAj k /dAu — Sji5 k i. Finally, we derive: 



(C5) 

(C6) 

(C7) 
(C8) 
(C9) 

(CIO) 

(en) 

(C12) 
(C13) 
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fdAEiA 



J = ^ l iEi t l /1 (A ll kSjm + A^jSkm) + £ivEi y mv(AvkSjl + A v jSkl), 

Im ' j k 



V dKlm / jk 

if I 7^ m, and: 
ZdAEiA 



= Ti^SkiAj^Ei^i + HvSjiEijvAvk, 



(C14) 



(C15) 



V dAu / jk 

otherwise. Writing the covariance matrix as V = V (0) +V (1) , the zeroth order reads V$° = {J2 l {x l -x {0) ) J (x l -x<- 0) )k/S l )/{T,d/S l 
The first order comes from the terms of Eq. CIO linear in Ei\ 

(i) 



V, 



En - 2E, u au.utA':°}E. 



(C16) 
(C17) 

where E = (E^/ /(EJ/ 'Si) and a jL „ v = {^{x, - x^)^ - x^E^ / S t ) / (^E,^ / Si). Note that V (1) = -E if there 
is no spread in the errors (i.e., Ei = constant). 



V [m = Elm — EfcjjO-m^ijjA^-Epi — TjkvOlk.mvA^l E vm 



APPENDIX D: THE PHOTOMETRIC ERROR MATRIX 

The photometric part of the error matrix Ei of Eq. |26| can be derived as follows. Paper IV shows that the errors in the 
half- luminosity radii Sri and average surface brightness Sui derived from the fits to the photometric profiles are correlated, 
with SFPi = Sri — aSui and a « 0.3, and calibrates the values of Sri and SFPi in terms of the quality parameter Q. 
Therefore, we consider the rotated coordinates (y, z) — T(log R e ,i, (5B e ,i)), where y = (logi? e ,i — a(<SB e ,i))/\/l + oi 2 and 
z = (alogi? e ,i + {SB e ^))/\/l + a 2 , and T the corresponding rotation matrix. In these coordinates the error matrix E yz is 
diagonal, with i? HZ ,n = <Jy and E yZ: 22 = c z . From the definition of y it follows a 2 = SFP 2 /(l + a 2 ). From logP e> i = 
(y + az)/\/l + a 2 and the fact that the errors on y and z are uncorrelated, we get Sr 2 = (a 2 + a 2 a 2 )/(l + a 2 ) and finally 
a 2 = ((1 + ct 2 )8r 2 - SFP 2 /(1 + ct 2 ))/a 2 . The matrix T T E yz T gives then: 

2.22 

Ei,n = " 2 = fr 2 , Dl 
1 + or 

P -P a (.-r 2 2,_ {l + * 2 )5r 2 -8FP 2 

Ei,is - E^i - (<x. ~ *») q(1 + q2) ■ (°2) 

a 2 ag + a 2 (a 2 - 1)<5FP 2 + (1 + a 2 )5r 2 
= 1 + Q 2 = a 2 (l + Q 2) ■ (D3) 



APPENDIX E: THE NORMALIZATION OF THE 3-DIM DISTRIBUTION FUNCTION 

The integral of Eq. ^ can be performed by defining the set of rotated coordinates (u>i , W2 , W3 ) = Q(x): 

x\ — bx:j 



Wl 



yr+fe 2 

W2 = X 2 , (El) 

bxi + x-s 



•U'3 



(E2) 

Writing M = Q T AiQ, integration over W3 leads to the following integral: 



oo poo 

T, 



f 3i = ^J- I I eyL V {-w 1 Ww/2)dw 1 dw2i (E3) 

" w\ c J W2c 

where the 2x2 matrix W is 

W<,j=Mi,.7 rr^> (E4) 



andw k = (FP C ut-FP)/Vl + a 2 , where FP cu t = 0.78 log PWt - 0.61 (Paper III, a = 0.3 « 6), FP = logi? e - <5j - a{SB e )) 
and W2c = logo-cut — logo-. Finally one gets: 

hi = L(h,k,p), (E5) 



where F(/i, fc, p) is the bivariate probability integral (Abramovitz & Stegun 1971), h — wicy/W2,2(l — p 2 ), k = W2cy/Wi,i(l — p 2 ) 
and p = — Wi,2/y / Wi,iW2,2. 
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